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Abstract: In the present work we derive a Central Limit Theorem for sequences 
of Hilbert-valued Piecewise Deterministic Markov process models and their 
global fluctuations around their deterministic limit identified by the Law of 
Large Numbers. We provide a version of the limiting fluctuations processes in 
the form of a distribution valued stochastic partial differential equation which 
Cd I can be the starting point for further theoretical and numerical analysis. We also 

present applications of our results to two examples of hybrid models of spatially 
extended excitable membranes: compartmental-type neuron models and neural 
fields models. These models are fundamental in neuroscience modelling both 
for theory and numerics. 

I _^ , Keywords: Piecewise Deterministic Markov Processes; infinite-dimensional 
vQ ' stochastic processes; law of large numbers; central limit theorem; global fine- 
ly^ , tuations; Langevin approximation; stochastic excitable media; neuronal mem- 
l' ' brane models; neural field models 

O ■ MSC 2010: 60B12; 60F05; 60J25; 92C20; 

m 

1 Introduction 

•i-H 

\^ I The present work studies the global fiuctuations of hybrid processes also called Piecewise Deter- 

$H ' ministic Markov Processes (PDMPs) on Hilbert spaces around a macroscopic limit. That is, we 

present a limit theorem for the rescaled fiuctuations of PDMPs which admit a deterministic limit 
in the sense of a weak law or large numbers. We prove that under the suitable convergence of the 
initial conditions the rescaled fluctuations converge to a generalised Ornstein-Uhlenbeck process 
given by the mild solution of an Hilbert-valued stochastic evolution equation. PDMPs generalise 
continuous time Markov chains and possess a wide range of applications. The interest in these 
processes results from the fact that they combine continuous deterministic evolution with random 
instantaneous events which entails their specification as hybrid processes. Additionally, PDMPs 
arise when non-homogenous continuous time Markov chains are turned into homogeneous Markov 
processes by the state-space extension [53]. Hence, our results also provide limit theorems for 
non-homogeneous Markov chains. The central limit theorem in the present study builds up on 
previous work [57| where we established a Law of Large Numbers for PDMPs and a limit theorem 
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for their internal fluctuations. The theorem in this study thus completes the 'trias' of limit the- 
orems as have been previously considered for time homogeneous Markov chains [211 [22] , PDMPs 
in finite dimensional Euclidean space [231 and reaction-diffusion models [SJ [T^ . 

In recent studies, PDMPs proved particularly useful for modelling in neuroscience and physiology 
under bottom-up ([3J [51 [IHl [13|) or multi-scale approaches ([121 [HI [21])- Indeed mathematical 
models for biological real-life processes are constructed on largely varying temporal and spatial 
scales. Depending on the perspective what seems to be a regular behaviour on one level of 
description arises as the emergent effect averaged over a large number of individual dynamics on 
a lower level at best described in a stochastic way. Examples in neuroscience arc travelling waves 
on excitable membranes which arise due to the gating dynamics of huge numbers of individual 
ion channels immersed in the membrane (cf. |18| ) or macroscopic models of brain activity where 
the activity is due to a huge network of interconnected individuals neurons (cf. |6]). However, 
when the macroscopic averaged pattern is observed then clearly not all stochasticity is lost. Even 
the higher level dynamics show random fluctuations albeit often typically small. Nevertheless, 
these small fluctuations can be of fundamental functional importance to systems characterised 
by complex non-linear responses which is the often the case in real-life systems. Hence, small 
fluctuations can cause tremendous changes in functionality, in the positive as in the negative way. 
For instance previously mentioned ion channel fluctuations can cause propagation failure of nerve 
signals and thus noise places limitations on the size of nerve fibres (cf. |14|). In many well known 
examples, fluctuations arc shown to be inversely proportional to the square root of the system 
size. In the present spatio-temporal framework however the difficulty is precisely to define what 
the system size is. This is an important point since determining the noise scaling constitutes 
one of the practically relevant applications of the kind of limit theorems we prove in the present 
study. Another application is the provision of a limit model of the fluctuations process that is 
supposed to be more tractable than the original process. This procedure is usually called diffusion 
or Langevin approximation. For the models in this study, this tractable model takes the form 
of a stochastic partial differential equation. However, as being distribution valued, it forbids a 
direct approach by standard approximation methods and new interpretation has to be provided 
first. We believe that the stochastic partial differential equation which we provide as a version of 
the limiting fluctuations processes can be the starting point for further theoretical and numerical 
analysis. This mathematically challenging and practically highly relevant point can be addressed 
on the basis of the present work. 

The remainder of the paper is organised as follows. In Sections [2] and \3\ we review the class of 
PDMPs used in this study as well as results we obtained previously which are relevant here. In 
particular, we present a law of large numbers for PDMPs providing a macroscopic 'average' model 
around which the fluctuations are studied. In Section [H we state the Central Limit Theorem 
which is the main result of this study and discuss its sufficient condition. The proof is presented 
in Section [5l In Section [6] we apply the limit theorem to PDMP models of excitable membranes 
and neural fields. The Appendix contains technical proofs of auxiliary results for the establishment 
of our main theorem. 



Remark 1.1 (Notation). We gather here notations regarding Hilbert spaces. We use * for the 
dual of an Hilbert space as well as the adjoint of a linear operator. Round brackets (•, •)/f denote 
inner product in the Hilbert space H. Angle brackets (•,-)h denote the duality pairing in H, 
i.e., the application of a linear, bounded functional in H* (first argument) to an element of H 
(second argument). Square brackets [•] are used to denote the elements where partial derivatives 
are evaluated at. Given two Hilbert spaces X, H the notation X ^^ H (continuous embedding) 
means that X d H and the embedding operator is continuous: \\u\\h < Cj|u||x for all u G A", for 
some finite constant C. An evolution triplet features X ^^ H ^^ X*, where X and H are Hilbert 
spaces; the inner product in H is identified with the duality pairing in X: for all G i7 and u € X, 
{(j),u)x — (0, "")//• These embeddings must be distinguished from the canonical identification of 



the Hilbert space X with its dual, therefore we use for this a distinguished notation: lx for the 
canonical embedding of an element on X into the dual X* and t^ for the Riesz Representation 
of an clement in X* by an element in X. Accordingly, (ti/(a;),y)x = {x,y)x for all y € X and 
L^ {4>) is the unique clement in X such that {(f>,y)x ~ (i^x {4>)Ty)x for all y ^ X. 

2 Spatio-Temporal Piecewise Deterministic Processes 

We give a brief definition of hybrid models / PDMPs relevant for the present study and refer to 
Uni HZl [2S] for a more detailed discussion. Let (fi, J^, (J"()t>o,P) be a filtered probability space 
satisfying the usual conditions, X and H denote separable Hilbert spaces forming an evolution 
triplet X C H C X* and K be an at most countable set. In this study all spaces and sets are 
equipped with their Borel-cr-fields and measurability always means Borel-measurability. Then a 
PDMP {Ut, 9()t>o is a cadlag strong Markov process taking values in H x K which is uniquely 
defined by the quadruple {A, B, A, /i) in the following sense: 

(i) The operators A : X x K ^f X* , which is linear in its A-argumcnt, and B : X x K ^f X* 
are such that the abstract evolution equations 

u = A{e)u + B{u,e) yeeK (2.1) 

arc well-posed in the weak sense for any initial condition uq G H, i.e., there exists a unique 
weak solution u to dm such that u € L^{{0,T),X) n H^{{0,T),X*) C C{[0,T],H) for 
all T > with u{0) = uq in H. We denote by 0((uo,^) the value of the unique solution 
corresponding to the parameter 9 at time t started in uq. Then the PDMP {Ut,Qt)t>o 
satisfies Ut = (jyt-TkiUrk) ior t G [Tk,Tk+i), where the random variables Tfc, fc G N, denote the 
jump-times of the PDMP with tq = 0. 

(ii) The measurable map A : H x K ^ M+ is locally path integrable along the solutions of (|2.ip . 
i.e, 

i-T 

A{(l),{u,9),9)ds < oo y{u,e)eHx K,T>0, 
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and defines the distribution of the jump times via 



'[rfc+i > 1 1 J-,,] = exp(- / A(</>, ([/,,, e,J, e,J ds) Vi > r^ . 



ft-Tk 

/o 
Hence, the process A{Ut, Qt) states the random instantaneous jump rate of the PDMP. 



(iii) Finally, the Markov kernel ji from H x K into K defines the conditional distribution of the 
post jump values, i.e., 

p[([/.,,e.jGA|([/.,_,e.,_j] ^ (5u^^_ x^((c/.,_,e.,_j,-))(A) 

for all Borel sets A m. H x K, with /i(([/T-j^_, 0t-^._J, {C/t-^-} x K) = 1. This implies that 
the components Ut and 0* possess continuous and piecewise constant paths almost surely, 
respectively. 

In this study we assume that all PDMPs are regular, i.e., Wnik^ao tu ~ oo a.s., and thus possess 
only finitely many jumps in any finite time interval almost surely. This assumption is particularly 
satisfied if the jump rate A is bounded. The domain of the extended generator of the process and 
its form for certain important functions is given in the next theorem, which is proven in [8l I25| . 



Theorem 2.1. (a) A hounded, measurable function f : H x K ^- R is in the domain of the 
extended generator of a PDMP if the mapping t !-> f{Ut, Ot) is absolutely continuous almost 
surely and the mapping (^, s,w) i— >■ f{Us-,£,) — /(C/s_(w), 0s-(ti;)) is integrable in the mean 
with respect to the random measure A{Us-, 8s-)/x((t/s_, Qs-), d^)ds. 



(b) Moreover, if in addition f is continuously Frechet- differentiate with respect to its first 
argument and such that the Riesz Representation /„ & H of the Frechet derivative satisfies 
fu{u,0) £ X for u G X and is a locally bounded composition operator in L^((0,r),X), then 
the extended generator Af is given by 

Af{u, 6) = {A{0) u + B{e, u), fu{u, 0))^ + K{u, 0) J (/(«, C) - f{u, 0)) m((w, ^), d^) ■ (2.2) 



An immediate consequence of Theorem 12.11 is Dynkin's formula for PDMPs. If the stochastic 
integral of / with respect to the associated fmidamcntal martingale measure is a martingale, then 
it holds that 

ni{Yt)m = f(Yo)+E\ f Af{Y,)ds Yo] . (2.3) 

'-Jo J 

Remark 2.1. Important for our study are the functions / of the following type 

f:HxK^C:iu,e)^ cxp{i (t^, u)h + i ($, ^(e))^) , (2.4) 

where E is a separable Hilbert space, i is the complex element, ip E X ^^ H* , z : K ^ E is a 
measurable map and ^ £ E* . We show in Appendix lA.ll that the function (|2.4p indeed satisfies 
the conditions of Theorem l2. 11 (b). 

2.1 A particular class of PDMPs 

Motivated by applications, we subsequently work with specific sequences of PDMPs that wc now 
describe. In applications, mainly the PDMP's continuous component U is the variable of interest 
which may be experimentally observed and one is usually only interested in this part of the 
information the piecewise constant component provides that is sufficient to define the dynamics 
of the continuous component. In the following we consider PDMPs indexed by n G N defined 
by the quadruples (A",i3", A",/i") on the probability spaces (fi", J"", (J"")t>o,P"), where we 
assume that there exists a separable Hilbert spacqj E and a sequence of measurable functions 
z"- : K"- ^ E, such that 

.4"(6i") = .4(z"(6i")), B''{u,9)^B{u,z''{9")) Vr e i^",n e N, 

for some operators A : E ^ L{X,X*) and B : X x E ^ X*. That is, the continuous dynamics 
depend on the piecewise constant component only via the functions z". Thus, actually, the limit 
theorems we derive for are limit theorems for the sequence (C/",2;"(0"))t>o (see [17]). 

Towards applications, we are particularly interested in excitable membranes and neural fields. The 
aim is to model on the one hand a spatially extended neuronal membrane with a finite number 
of channels immersed and on the other hand a finite populations of neurons in a brain region. 
Indeed this corresponds to biologically realistic situations and to experimental conditions as well 
since in an imagery treatment only finite number of neurons can be observed. For these two 
cases, using the spatial character of the models, we consider partitions of a full spatial domain 
D C M."^ into subdomains Z?*"''", k = 1, . . . ,p{n). In the case of an excitable membrane equipped 
with ion channels this is a partition of the membrane patch. Whereas in the case of neural 
fields it is a partition of the brain area of interest. Then the piecewise constant component 
9^ = (6*^'", . . . ,6'P("^'") G K^ C N^(") counts the number of open ion channels (or the number 
of active neurons) in the subdomains D*^'". The total number of channels (or neurons) in the 
subdomain D*^'" is denoted by l{k,n) < oo. In these two cases the coordinate functions z" are 
given by 

P(") nk.n 

^''in-T.j^^^n.^^L^-iD), (2.5) 



^In the study 1271 the functions z", called coordinate functions, arc actually vector- valued with components in 
a Hilbert space and then E corresponds to the direct product of the spaces which is a Hilbert space itself, hence 
no generality is lost presently. 



which gives at each point in space the fraction of open ion channels (active neurons) in the 
respective domain. Due to the spatially piecewise constant nature of the coordinate functions in 
we call these PDMPs Compartmental Models. 

2.2 Quadratic variation associated with PDMPs 

Before proceeding further, we now discuss a fundamental fact regarding the aim of deriving limit 
theorems in our framework; the decomposition 

z"(en = z"(e^') + / (yiz")(?7f , e^dt + Af^ (2.6) 

of the piecewise constant component of the PDMP. The integral in the right hand side of (|2.6p exists 
in the sense of Bochner and the process M" is a cadlag, square-integrable S-valued martingale 
under Assumptions |X] displayed in the next section. As (|2.6p can be considered a stochastic 
evolution equation driven by the martingale M", one can think of these martingales containing 
the inherent stochasticity in the process. To prove a central limit theorem one has to investigate 
more closely the structure of the martingale M" in (j2.6p and establish a non-trivial limit under 
suitable rescaling. 

The quadratic variation G" G L{E,E*) associated with the PDMPs (?7",8") is the mapping G" 
defined via a bilinear form on £ given by 

($,vi/)K^(G"(f/r,er)$,*)£ ^ 

A"(t/r, er) j ^ ((z"(0 - ^"(er), $)£) ((^"(C) - ^"(er), *)£) ^l-{(ul\ e?), dc). (2.7) 

Then it holds that the operators G^{U",z'^) are symmetric, positive and, anticipating the condi- 
tion (jC3.ip below, of trace class. Furthermore, for the martingale part in (|2.6p it holds 

E"\\Ml'\\l. = E" / TrG"(C/f,endi (2.8) 



E" 





t 



A"(c/r,er) / \\z"iO - ^"(eDi. /."((c/;\er),dO di 

Jif" 



For a more detailed discussion of quadratic variations associated with PDMPs and proofs of the 
cited results we refer to [27l Section 3]. 

3 Law of Large Numbers 

In this section we review the Law of Large Numbers obtained in |27| which establishes the con- 
vergence of the PDMP sequence to a deterministic trajectory. Then we present a simple corollary 
relevant for the present study. 

In our framework the deterministic limit solves an evolution problem as follows. Let F : Xx E ^>- E 
be an operator such that the deterministic abstract evolution system 

ii = A{p)u + B{u,p), 

(3.1) 
P = F{u,p) 

is well-posed in the following sense: for suitable initial condition (mo,Po) & H x E the system 
p.ip possesses a unique global weak solution such that for all T > the component u is in 
H^{{0,T),X*)nL^{{0,T),X) C G([0,r],iJ) andpisin H^{{0,T),E) C C{[0,T],E). We now 
describe the set of assumptions (Assumptions |^ under which we proved the Law of Large Num- 
bers. 



Assumptions A. We assume the following conditions arc satisfied: 



(Al) The operators A, B and F satisfy almost surely Lipschitz-type conditions uniformly in n 
along the deterministic and stochastic solutions, i.e., for all T > there exist constants 
ii, L2, independent of n, such that almost surely 



(A(z"(en) C/r - A{p{t)) u{t), [/," - u{t))x + (i?(C/", z"(e")) - B{u,p), [/" - u)x dt 



< ii / \\ur-u{t)\\% + \\z"{e-)^p{t)\\ldt. 
'0 



and 



\\F{u,p) - F{U\z^&^))U^io.r},E) < L2\\{u,p) - {U",z^e"))U^^o,T),H.E) 

Note that the constants Li, L2 niay depend on T and {u,p). 
(A2) The initial conditions converge in probability, i.e., for all e > 



lim 



\\K-uorH + \\z^{e^)-po\\l>e 



= 0. 



(A3) The fluid limit assumption holds in probability, i.e., for all T, e > 



lim I 

n— ^cxD 



\[A"{-,z-i-))E]{u:,e:)^F{u:,z"{e:))\\ids>e 



(A4) The second moments of the jump sizes vanish in the limit, i.e., for all T > 



lim 



A"([/:,e^) / ||z"(0-^"(e:)|||/."(([/r,e,"),dO 



ds = 0. 



Theorem 3.1 (Law of Large Numbers H?! Thm. 4.1]). Under the Assumptions [^ the se- 
quence (U", z"{&"))t>o converges uniformly on compacts in probability (u.c.p.) to the solution 
{u{t),p{t))t>o of dniH), that is, for all T, e > 



lim 



sup*e[o,T] iic/r-«(i)iiH + ii^"(0r)-pWi 



> e 



= 0. 



(3.2) 



Heuristically, the proof of Theorem 13.11 (law of large numbers) is done by inserting into the de- 
composition p.6p the limits for the initial condition and the integral given in assumptions (i^ 
and (A[3]) and noting that the martingale M" in (|2.6p converges to zero due to assumption (A|3]). 
The resulting equation is just the equation for p in (|3.ip . 

Remark 3.1. The u.c.p. convergence p.2p immediately extends to any normed spaces the Hilbert 
spaces H and E are continuously embedded in. Moreover u.c.p. convergence implies that the 
processes also converge weakly on the space of cadlag functions D(R+, H x E). 

Conversely to the extension of the law of large numbers to spaces with weaker norms, it can also 
- albeit not immediately - be extended to spaces of higher spatial and/or temporal regularity, 
provided that the continuous component of the PDMP sequence and the deterministic limit are 
sufficiently regular. The extension of the results works along the lines of standard estimation 
techniques from the analysis of (linear) evolution equation (see, e.g., [IS])- (For a specific model, 
a more direct method was used in [3].) A first extension, which we employ later in this study, is 
stated in the following corollary. 



Corollary 3.1. Let T > be fixed and assume that the operator A is independent of p and satisfies 
the energy estimate 

li\\u\\\<-{Au,u)x+l2\\u\\]j yueX 

for constants 71 > and 72 > and that the mapping B : H x E ^ X* is continuous. Then for 
all T, e > 0, 

lim P"[|l[/"-7.|U2((o.T),x)>e] =0. 
Proof. By dcfinitiorG of a weak solution, the following holds almost surely for almost all t G [0, T], 

({/,"- u(i),^7^-^i(t))x 
= {AUr- Auit), Ur - u{t))x + {B{Ur, z"(en) - Biuit),p{t)), Ur - u{t))x . 

The energy estimate yields 

7irr - u{t)\\\ < 72rr - ^mi + iis(c/r,^"(er)) - B{u{t),pm\x4ui^ - umx 

and due to Young's inequality we get for all e > 



(71 - e) \\u- - umi < i2\\ur - u{t)\\i + \ i!s(t/r, z-ien) - B{u{t),p{t))\\ 



Choosing e small enough such that the left hand side of this inequality is positive and integrating, 
we find that 

W - Al^mT).x) <K\\u-^ Almo,T),H) + K ||i?(f/", z"(e")) - i?(«,p)|lL((o,T),x*)> 

for some eonstant K. Here the right hand side is a random variable that converges to zero in 
probability due to Theorem l3. II and the Continuous Mapping Theorem. The corollary now follows 
immediately. D 

From Theorem l3.1l we see that all stochasticity of the PDMP sequence is lost in the limit since this 
theorem provides a connection of the sequence to a deterministic, macroscopic equation. At this 
point two types of fluctuations come to mind: one can either consider only the internal fluctuations 
of the piecewise constant component z"(0") or the global fluctuations PDMP sequence. In p7] 
the internal fluctuations were characterised via a central limit theorem for the martingale part of 
the PDMPs, i.e., the process M" in (|2.6p . In this study we consider the global fluctuations of the 
PDMP sequence around its deterministic limit. 

4 Central Limit Theorem for global fluctuations 

In this paper we investigate the fluctuations of a PDMP sequence around its deterministic limit 
given by the law of large numbers of Theorem 13.11 Hence we are naturally concerned with the 
differences U^ — u and z"(6")— p which we have to rescale in order to obtain a non-trivial limit. We 
thus introduce a real, positive sequence a„, n G N, necessarily diverging, and make the stochastic 
process given by 

z," := v^(c/r - «(i),^"(er) -pW), t > 0, (4.1) 

the central object of interest in our study. Such a rescaled fluctuations process is classical in the 
study of central limit theorems. The main result we presently establish is that the sequence of 
processes (|4.ip admits a converging subsequence with limit supported on the space of continuous 
functions and with a further strengthening the limit is a unique diffusion process, see Theorem 



^It holds in between jumps for almost all t and as the process is regular, there are almost surely only finitely 
many jumps. Thus for each u) the possible set, it does not hold, is a finite union of Lebesgue null sets, hence a null 
set. 



It is important to note that for central limit theorems in infinite-dimension it is in general not 
possible to prove the convergence in the norm of the Hilbert space H x E for which the Law of 
Large Numbers holds, nor docs the limit exist in this space. However, we can obtain convergence 
in weaker norms: in the literature limits for spatio-temporal models are often considered in the 
Schwartz space of distributions. We consider the limits in duals of subspaces of H and E. A 
natural choice for the subspace of H is the space X in the evolution triplet. The convergence will 
take place in its dual X*. In order to keep the presentation simple we restrict ourselves to this 
choice. For the limit of the coordinate functions, wc introduce a Hilbert space £ ^^ E usually 
obtained by the completion with respect to a different inner product defined on a subspace of E. 
Clearly, together with the duals these spaces form an evolution triplet with E C £* and the norms 
satisfy ||z|J£» < CHzJI^; for all z G E for some constant C. In the application to neuroscience 
models E := L^{D) and £* is given by the dual to some Sobolev space H^{D). For a more 
detailed discussion see Remark 16.21 When establishing tightness of the sequence (|4.ip we assume 
for technical matters that there exists a third Hilbert space V satisfying the chain of continuous 
embeddings 

£'^V'^E'^V*'^£*. 

In applications these embeddings of spaces may collapse, i.e., some or all the spaces coincide. This 
is indeed the case with finite-dimensional Hilbert spaces where all the above spaces coincide with 
an Euclidean space. This property renders limit theorems in finite dimensions much simpler. 

4.1 The main result on global fluctuations 

We now proceed to our main result. We will establish the Central Limit Theorem under the 
following set of assumptions that we call Assumptions C. 

Assumptions C. (CI) The operators A, B and F satisfy Lipschitz-type conditions, uniformly 
in n, along the stochastic and deterministic solutions: for all T > 

(A(z"(en) C/r - A{p{t)) u{t), [/," - u{t))x + {B{U", z"(e")) - B{u,p), C/" - u)x dt 

< lJ \\u--u{t)\\i + \w\Q'i)~p{mi'dt, 

Jo 

\\F{u,p) - F{U", z"(e"))|U2((o,T),y.) < L,\\{u,p) - ((7", z"(e"))||L^((o,T),HxV.), 

\\F{u,p) - F{u-, z-{en)\\Lmo.T).£') < L,\\ {u,p) - ([/", z-ien)\\Lmo.T)Mx£') 

almost surely for some constants L3, L4 and L^ which do not depend on n but may depend 
on T or {u,p). 

(C2) The sequence of initial conditions satisfies 

supV^E"||t/o"-7.olU + sup V^E"||z"(e^) -Polk- <«). 

neN nefi 

(C3) For all T > 0, 



supa„ 

riGN 



A"(t/:\e:)/ ||z"(0-z"(e;')lk./^"(([/r,e';),de)ds < 00, (C3.1) 



and 



riGN JO 



T 



supa„ / E"||[^"(z7(.), •)y](t/r,er)-F([/:,z"(e:))f^.d,s < 00. (C3.2) 



(C4) There exists an orthoiiornial basis {ifk)k€fi oi X x £ such that for all S > Q and almost all 

lim supP"[V |(Zr,^fc)xxd' > ^1 = 0. 

(C5) For all $ e £ and all T > 0, 

lim / E"|(GKs),p(s))$,$),-a„(G"(C/;,e:)$,$)Jds = 0, (C5.1) 

where G'(m(s),p(s)) G L{£,£*), s > 0, is an integrable family of symmetric, positive, trace 
class operators, and 

lim V^/ E"|[^(z"(.),$)£]([/r,e';)-(F(C/;,z"(e:)),$)Jds = 0. (C5.2) 

ri— i-oo Jq 

(C6) For all $ e f and all T, c> 0, 



lim E" / A"([/;,e^) / M"((C^r^0"),de)ds =0 (C6.1) 

and the jump heights of the rescaled martingales are almost surely uniformly bounded, i.e., 

supsup V^||z"(e't') - z"(e't'_)||£. < oo a.s. (C6.2) 

neN t>n 

(C7) The operators A : X x E -)■ X* , B : X x E -)■ X* a.nd F : X x E -)' E arc partially 
continuously Gateaux differentiable along the solution {u,p) of the deterministic system and 
for almost all s > 0, the Gateaux derivatives evaluated at {u{s),p{s)) possess a bounded, 
Lipschitz continuous extension to X* , respectively £"*, which we denote by ~. For instance, 
if Bu[u,p\ £ L{X,X*) denotes the partial Gateaux derivative of u i— >■ B{u,p) evaluated at 
{u,p) E X X E, then Bu[u,p] is its extension to e L{X*,X*) such that 

Bu[u^p\h ~ Bu[u,p]h yh£X. 

(C8) For some p > 2 and for all T > 0, 

sup E"i|C/,"-u(i)||P.+ sup E"||z"(ej')-p(t)||| < oo. 

nGN, t<T nen,t<T 

We now state the main result of this study. Remember that 

zr := ^/^([/r - w(t), z"(er) - pw) , i > o . (4.2) 

Theorem 4.1 (Central Limit Theorem for the global fluctuations). We assume that the Assump- 
tions |21 awrf O are satisfied. Then the laws of the processes (Z")t>o are tight in iD(R+,X* x £*) 
and any accumulation point P* is supported on the space of continuous functions, i.e., 

P*(C(M+,X* x£*)) =1. (4.3) 

Furthermore, the limit is unique, i.e., the processes {Z")t>o converge in distribution to an X* x£* - 
valued process with characteristic function tp : M^- x X x £ — > M satisfying the partial differential 
equation 

^ = -\{G{u{s),p{s))^, $)^^+(||,I*(p(s))^ + B:Ks),p(s)]^. + KKs),p(s)]$)^ 

+ (^^^A;[u{s),p{s)]i^ + B;[u{s),p{s)] i, + F;\u{s\p{s)\ $)^ (4.4) 

where Lp and its partial derivatives are evaluated at {t,ip,'^), with initial condition (^(0,V'j^) = 
lim„_j.oo E"e*''^o ''■''''*" ■^^'^ . Moreover, if the initial conditions Zq converge to a Gaussian random 
variable on X* x £* , then the limit is a Gaussian process. 



Remark 4.1. It is easy to see that Theorem 14.11 is an extension to previous resuhs considering 
PDMPs on Euchdean spaces. In finite dimensions hnear operators arc matrices with the adjoint 
operator given by the matrix transpose and thus equation (j4.4p reduces to the equation reported 
in [53]. Furthermore, reducing to piecewisc constant finite dimensional processes, i.e., A = B = 0, 
we obtain the equation as in Kurtz |22| . 

We stated that in previous work we presented a limit for the internal fluctuations of the PDMPs 
(cf. [27l Thm. 5.1]). The present Assumptions ICl are stronger than the assumptions needed for 
that limit theorem (if we assume that the embedding of £ g 1/ is of Hilbcrt-Schmidt type as in 
applications is always the case). Hence under the Assumptions [Cj it follows that the f* -valued 
martingales (Y^a^M")t>o converge weakly to an f *~valued centred, continuous Gaussian process 
with independent increments characterised by the covariance operator 

C(i) = / G(M(s),p(s))oi-ids Vt>0. (4.5) 

Finally, in [5S1 [57] a second version of the martingale central limit theorem with slightly different 
conditions was considered, i.e., a marginally stronger condition instead of (|C6.1[) and a marginally 
weaker condition instead of (|C6.2[) . We note that Theorem 14.11 remains valid if the Assumptions 
[Cl are changed accordingly. 

It is important to note that if the limit is a Gaussian process, it may be identified with the mild 
solution of a linear stochastic partial differential equation with additive noise, i.e., a generalised 
infinite dimensional Ornstein-Uhlenbeck process. 



Corollary 4.1. Assume that Theorem 14.11 holds and the initial conditions Zq converge to a 
Gaussian limit (Uq,Pq) e X* x £* . Then the limiting process is a version of the mild solution to 
the linear, non- autonomous stochastic evolution equation 



dUt 
dPt 



{A{p{t)) + Bu[u{t),p{t)])Ut + Fu[u{t),p{t)]Pt 



dt 



{Ap[p{t)]+Bp[u{t),p{t)])Ut + Fp[u{t),p{t)]Pt dt + g{t)dWt 



rQ 



(4.6) 



with initial condition {Uq, Pq), where W'^ is a Wiener process on a Hilbert space H' with covariance 
Q&L{H') andg{t) e L{H',r) such that g{t)Qg*{t) = G{u{t),p{t)) o q^ . 

Remark 4.2. It is always possible to find operators g in (|4.6p . For instance, choose U ~ £* 
and Q ~ I being the identity. Then g{t) is the unique square root of the trace class operator 
G{u{t),pit))oi-\ 

4.2 Discussion of the assumptions 

Before proceeding to the proof of Theorem 14. 11 we discuss the set of Assumptions G that we have 
chosen to rely on. 

4.2.1 Proving tightness under Assumptions [C] 

In the sequel we prove convergence in distribution of a sequence of cadlag processes. This is a well 
established question in Probability Theory for which sufficient conditions as well as the plan of 
the proof are well identified: one first proves tightness of the given sequence and then shows that 
the set of accumulation points is a singleton using some characterisation. Tightness of probability 
measures on the space of cadlag functions taking values in X* x £* is equivalent to tightness of the 
sequences of marginals (in a dense subset of M+ ) and convergence in probability for the mode of 
continuity. These properties are classically obtained by checking the conditions (TlT])"(TSl) below. 
The first property is implied by (Tjl]) and (T^ (which are also necessary) , the second one is implied 
by condition (T^ due to theorem [H Thm. 8.6, Ghap. 3]. 
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(Tl) There exists an increasing sequence of finite dimensional subspaces 5'fc, A: £ N of X* x £* 
such that hmfc_>.oo T^kX = x for all x G X* x £*, where tt^ denotes the projection onto Sk, 
and for every e, 5 > and all t out of a dense subset of M+, there exists Sm out of (S'fc)fcgN 
satisfying 

supP"[inf^eS„J|Zr-2/|U.x£. ><?] < e ■ 



ri6F 



(T2) For all e > and all t out of a dense subset of R+, there exists a (5 > such that 



sup 

raGN 



Tin r II ryn 



Zl'Wx'xe- >6]<e. 



(T3) For every T > 0, all At G (0, 1) all n e N, there exist non-negative random variables J" (At) 
such that for all < i < T and < s < At 



E" 






IX'Xf 



^r 



< 



E"[^"(At)| J-j" 



and 



lim limsupE"[J"(At)] =0. 

At->0 ,i_j.oo 



(4.7) 



(4.8) 



However in the present paper we found it more natural to work with the set of Assumptions C. 
We now show that Assumptions C imply (TlT])-(T|3]): 

To this end we first note that if S'„ are the spaces spanned by an increasing set of elements of the 
orthonormal basis obtained from the Riesz Representation of the orthonormal basis {^Pk)k£'H in 
X X £ then, assumption (CH]) is just a reformulation of condition (lH])- To proceed we need the 
following lemma, which is also of later use. 

Lemma 4.1. Assume that (C[T]) holds. Then there exist constants Ki, A'2 independent of n ^H 
and t < T such that for all n €z N and t < T, 



a„E" 



< 



\\u-^u{t)\\l + \\z-{e-)-pml,_ 

(a„ E"||C/o"-uo||ff+a„IE"|k"(ejf)- Polly. +a„E" sup ||M;||2,.) i^^e^^* (4.9) 
^ se[o,t] ^ 

+ [a^E^ J)\A^{{;z-{-))e){U:\Q:) - F{U:,z^Q:))\\l,ds)tK,e''^\ 
Assume moreover that (C^ and {C^ also hold. Then the family of random variables 

{V^\m " ^{t)\\H + V^||^"(en -P(t)lks t e [0,T],ne n} (4.10) 

is uniformly integrable for every T > 0. 

Wc postpone the proof of Lemma 14.11 to the next section and prove now that (Tj^)) and (T[3]) 
hold under Assumptions C. We see that Lemma 14.11 combined with the continuous cmbcddings 
H '^ X* and V* ^^ £* , and the inequality y^o^+l?) < \a\ + \b\, implies that also the set 
{||Z"|jx*xf*; t G [0,r],n e N} is uniformly integrable. Thus there exists an increasing, non- 
negative, convex function g : [0, 00) — > M_|- such that lim2;_j.oo g{x)/x = c» and 



sup E"[5(||ZriU.x£.)] <oo, 

tG[0,T],nGN 

[T^ Appendix, Prop. 2.1]. Now, using Chcbyshcv's inequality we obtain 

1 



(4.11) 



'[\\zn\x'xe'>S]< 



9iS) 



E"[5(ll^, 



t WX-xS 



)] 
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Taking the supremum over all rt G N in the right hand side of this inequality and choosing S 
sufRciently large yields (TjS)). 

Moreover, using the norm oi H x £* yields 

M II ? 

= a„ \\ui^,-u(t + s) - t/r + ^iWIII/ + «„ ||2"(eivj -p(t + s) - ;^"(er) -p(t)|l|-. 

Next the estimation techniques used to prove [Ul Thm. 4.1] yield the inequalities 

\\\J'^^,~u(t + s)-Ul^ + um\ (4.12) 



< 4L., 



t+s 



\\U^ _ u{r) - U- + u{t)\\l + ||z"(e';) -p(r) - ^"(9?) - p(f)||2.. dr 



+ A.sL,[\\U--uit)\\l + \\z-ie-)^pmv'), 



and 



||z"(er+j-p(t + s)-z"(er)-pWil|. 



(4.13) 



< 6U 



t+s 



\\U- - u{r) U^ + u{t)\\l + \\z-{e-)-p{r) - ^"(6^) - p(t)li^ dr 



+ 3 WMlXs - Mill' + 35i5(||C/r - uml + |l^"(en -P(t)ll^) 



3s 



t+s 



|^"((z"(.),->£)(c/;\e;')-F(z"(e;'),c/;)g.dr, 



where L3 and L5 are the Lipschitz constants in assumption (CH)). Remember that G" is the 
quadratic variation defined in (|2.7p . Then, for fixed i > 0, the process 



s^iiA/,%-ii/rii|. 



t+s 



TrG"([/;',e;?)dr 



is a martingale with respect to {J-'l]^^g)s>o satisfying 



E" [111^,"+, - M^l, I JT] = E" / Tr G"(C/;\ ep) dr 



t+s 



TV 



This can be proven using the martingale results developed in [271 Section 3] . Now, a combination 
of the estimates (|4.12p and (|4.13p . employing the continuous embeddings £ ^^ V and the condi- 
tional Fubini Theorem, see, e.g., [31 Thm. 1.1.7], for interchanging conditional expectation and 
integration, yields 



E^ r ^ 1 1 1 ^7^ rr^ 1 1 2 I T-nl 

]."« W^t+s- ^t ll^xf*!-^* J 



t+s 



< i^i / E"[a-il|Z;'-Z«||l,xV.|-^r]dr + E" 



t+s 



TrG"(C/;",e;!)dr 



^r 



+ ,s A-i(E"||C/r - «(i)lll, +E"l|z"(er) -p(i)lPv-) 



s JsTiE" 



t+s 



|^"((z"(.), •)£)((7,",e;?) - ^(z"(e;'), [/;)||^. dr j-f 
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for a suitable finite constant Ki independent of n G N. We then apply Gronwall's lemma to this 
estimate and multiply the resulting inequality by a„. Thus, as s < At, wc obtain the estimate 

E"[iiz[v.-^riiffx£.|-^r] 

rt+At 

TrG"(t/;?,e;!)dr 






^K,At 



< an K 

+ At Ki an E" 



,^iAt 



V"((^"(-), •>£)([/;, e;?)-f(z"(e;?),c/;)||J.dr 



-pn 



^K,At 



Let us define the random variable S"{At) in (|4.7p as being the right hand side of this inequality. 
This choice is relevant since the left hand side is an upper bound to E"[||Z("[_g — ^"Hx-xf |-^"] 
due to the continuous embeddings H ^^ X* . Then this random variable satisfies 

pt+At 

limsupE"[,5"(Ai)] < i^2At + /v2 InnsupE" a„ / TrG'"([/,", 9^?) dr 

l-t+At 

< K2At + K2 TTGiuis),pis))ds (4.f4) 

for a suitable constant K2 independent of n. Here the first estimate is due to condition (CJH]) 
and Lemma (|4.ip . The second one follows from the convergence result in Proposition I A. II in 
the appendix. Since the upper bound (|4.14p converges to zero when At — )> 0, condition (TjSl) is 
satisfied. 

4.2.2 Norms and spaces 

In our assumptions wc have to choose the right norms and spaces carefully. In applications 
the choices for the correct spaces X* and £* are usually determined by the need to establish 
condition (CS|). If the embeddings X ^^ H and £ ^^ V are of Hilbcrt-Schmidt type then (C2]) 
is satisfied under the conditions (CH]) - (CH]), see Remark [CTI Thus, in applications one usually 
chooses first the space V* small enough such that conditions (C[T])-(C[3]) hold. Then the spaces 
X and £ are chosen large enough such that the embeddings X ^^ H and £ ^^ V are of Hilbert- 
Schmidt type. Choosing largest possible spaces V,X,£ one obtains the best possible regularity 
of the limiting process. Secondly, condition (jC3.1|) implies a condition akin to (i^ and thus 
the processes (^a„ M")t>o are squarc-intcgrable V*~ and £'*-valued cadlag martingales under 
this assumption. Moreover the assumption {C^ implies the same statements for the space £* . 
Concerning assumption (CO, we note that the required extensions exist if the embeddings are 
dense which is usually the case in applications. Also, in applications (C^ is cither a consequence 
of the Law of Large Numbers or follows directly from the path properties of the PDMPs. 

5 Proof of the Theorem 14.11 (Central Limit Theorem) 

We split the proof of Theorem 14.11 into several parts. In the first part. Section 15.11 we estab- 
lish tightness of the sequence of laws of Z" in (|4.ip and thus existence of a weakly converging 
subsequence. The short second part. Section 15.21 shows that any weak limit of a converging 
subsequence is a continuous stochastic process. The central part of the proof is the third one, 
Section 15.31 Therein we first prove that the characteristic function of the limit satisfies the par- 
tial differential equation (|4.4p . To this end we have adapted a proof from [551 US] to the present 
infinite-dimensional setting. The main difficulties here are technical in nature by identification 
and careful treatment of corresponding infinite dimensional concepts that substantially simplify 
in Euclidean space, that is, the use of general functional analytic operator theory instead of linear 
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algebra for matrices and vectors is necessary. We show that equation (j4.4p possesses a solution 
and the solution characterises a diffusion process. 

5.1 Tightness 

Wc now prove Lemma 14.11 Remember that in the preceding section this lemma was the key to 
prove that conditions (TIT|)^(TE1)j which ensure tightness, holds true under Assumptions C. Hence, 
tightness will follow. 

Proof of Lemma \4.1\ Due to the de la Vallee-Poussin Theorem [T2J Appendix, Prop. 2.2] it is 
sufHcient to prove that the family (|4.10p possesses uniformly bounded second moments. Using the 
first two Lipschitz conditions in (Cjl]) and Gronwall's inequality yields the almost sure inequality 



\U-^u{t)rH + \\z^{e-)-p{t)\\l,< (\\U--uorH + \\z-{e-)-porv^+ sup ||ACI|^.)i^ie^^ 

^ se[o,t] ^ 

+ (/'V"((-:^"(-)>v)(f/r,e^)-F(z"(e^),c/:)||J,.ds)i^ie^^* 





where the constants Ki, K2 are independent of n G N. Applying Cauchy-Schwarz inequality to 
the last term, taking expectations and multiplying by a„ yields (j4.9[) . Let us show now that each 
expectation in the right hand side of (j4.9p is bounded over all n <E N and t < T: Indeed the 
terms containing the initial conditions are bounded due to assumption (C^- The boundedness 
of the martingale term holds due to Doob's inequality and condition (|C3.ip as ||M"||y. is a 
submartingalc. The last term is bounded by condition (jC3.2[) . This proves the lemma. D 

The results so far have the following consequence which will be useful in the remainder of the 
proof. As we now have established tightness there exists a weak limit Zt to a subsequence of 
the sequence Z" in X* x £* . Without loss of generality we identify in the following Z" with a 
converging subsequence. Then, the Continuous Mapping Theorem and the established uniform 
integrability of Z" imply, due to [121 Appendix, Prop. 2.3], the convergence of expectations 

lim E"|jZ,"|U.x£. = ElJZtlU.xf. <oo Vt>0. (5.1) 



Remark 5.1. Before closing this section let us mention another consequence of Lemma [4~T1 Let 

fk = {fk'' "^f )' k E N, be an orthonormal basis oi X x £. The Markov inequality yields 

k>7n 

and employing the properties of the evolution triplet we further estimate 



k>m 

^ 1^" [E,>,„(ii^^ii?f + W'fi'rv) «" {\\ur - uimi + «« ik"(er) - pmi. 

Then the additional assumption that the embeddings X ^^ H and £ ^^ V are of Hilbert-Schmidt 
type implies that the sum J2k>m{\\'^k\\'H + llv'f llv) is finite and vanishes for m — > 00. Thus 

k~>m 

Lemma HTTI ensures that the supremum is bounded. Hence in this case (CS]) follows from (C[T|) - 

(cH). 
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5.2 Support of the limit (|i3| 

Let us show that any weak Hmit to a subsequence of the sequence Z" is a continuous process. This 
result is an immediate consequence of the Martingale Central Limit Theorem [571 Thm. 5.1] and 
its method of proof. Therein it was shown that the rescaled martingales converge to a continuous 
process and for this result it is sufhcient that 

lim E y a„||AA/t"|||. =E V ||AA/t|||. =0 (5.2) 

n— >-oo ^ — ^ ^ — ^ 

se(0,i] s£{Q,t] 

holds for alH > 0, where AfM" = Af" — M"_ and M denotes the weak limit of the martingales 
.y/a^Af ". The property (|5.2p was established in [57] under assumptions implied by the Assumptions 
[Cl An analogous result immediately follows in the present case for the processes Z". Note that 

||AtZ"||i-.^£. = aJMU" - u)\\j,, + an\\At{z"{Q") ~ p)\\l ^ + a„||AtA/"|||. . 

Thus 



lim E y \\AtZ 

s6(0,t] 



n||2 

A'*x£* 



follows immediately from (|5.2p . Then we infer from [T^l Chap. 3, Thm. 1.2] that the limit possesses 
almost surely continuous paths. 

5.3 Characterisation of the limit 

In the course of the proof we have so far obtained that there exists a weak limit to a subsequence 
of the processes Z", which is an a.s. continuous process on X* x £* . In this final part we now 
characterise this limit as a particular diffusion process and show that it is unique. We use a 
characterisation by convergence of the finite-dimensional distributions. Let us introduce some 
notation. We denote by {Zt)t>o ^ version of the limiting process defined on some probability space 
(17, T, P) and E denotes the expectation with respect to the measure P. For given ip £ X C H 
and 4> S f we denote by h the bounded, continuous mapping 

h:Hx£*^C:{x,y)^ cxp(i(x, V)x + i{y, $>£) (5.3) 

for which we also use the shorter notation h{z) = h{x,y) for z = {x,y). Moreover, due to the 
chain rule h is partially Frechet differentiable with derivatives 

dh dh 

— ix,y)^ihix,y)^ e LiH,C) = H, —{x,y) = ih{x,y)^ S L{8* X) ^ 8 . 

That is, the derivatives lie in the corresponding dual spaces, i.e., H* = H and £, respectively, 
considered as complex Hilbert spaces. Conversely, h can be considered as a function of V' and $ 
for given x and y in which case we obtain the partial derivatives 

^ J oil 

— {x,y)^ih{x,y)x e LiX,C)^X*, —{x,y)^ih{x,y)y e L{E,C)=£*- 

By the use of this notation the characteristic function ip"" of the random variable Z" e X* x 8* 
satisfies 

^"(i,7A,$) = E"/i(zr) = E'hi^^^iuit - M(t)), v^(2"(er) -p(i))) . 

Note here that t/" — u{t) is considered a random variable in X* , however, it actually takes values 
in H which is embedded in X* in the sense of evolution triplets. We show in Appendix (jA.l[) that 
h defined in (|5.3p is in the domain of the extended generator of the PDMP ([/", 0", tt(i),_p(i))t>olj 



^It is clear that the process ([/", G", M(t), p(i))t>o is a PDMP as the components (t/", 0")t>o form a PDMP and 
also (^(t), p{t))t>o are a (degenerate) continuous Markov process as being the solution of a deterministic abstract 
evolution equation. 
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Thus Dynkin's formula (j2.3[) yields |f| 



E-[{A{z-{e:))U: + B{U:,z-{e:))-A{p{s))uis)-B{u{s),p{s)),i;)^zV^,h{Z:) 
-{F{p{s),u{s)),'^)^t^W,h{Z:) 
+ A"{U:, &:) f (/^(^/^(C/^ - u{s)),^^{z-{0 - p{s))) 

- h{y^{u: - uis)), v^(z"(0:) - Pis)))) M"((f/r, er), d^)] ds 

= : / E"[ifi(.s)+i72(s)]ds, 

where the terms Hi and H2 denote the continuous and jump part of the generator, respectively. 
Introducing the function 

K{u) = e''"- 1-m+iu^ 

for which K{u) = 0{u'^) holds, we obtain 

H,{s) = A-{U:\e^)h{Z^) J [exp(v^(z"(C)-^"(e^),*>£)-l]A*"((t/r,0"),dO 

= A"([/;, &:)h{z:) f f*v^(z"(0 - z"(e»), $)£ 



-ffn 



ttn 



Next note that (see [27l Sec. 3] for details), 

A"iu^,&:)J ((z"(o,<i>)£-(z"(er),<i>)£)A*"((c^r>0")>dO = [^(z"(-),<i')£](f/r,e:), 

where y^(z"(-), $)£ denotes the application of the generator of the PDMP to the map ([/, 6, u,p) 1— >■ 

Therefore we obtain that (p'"^{t, i/.', $) — (yc"(0, 4', ^) is the sum of five integral terms that we denote 
by (i) to (v). Next we make these terms explicit and we discuss their respective limits for 71 — > 00. 
We start with (i), (iv) and (v) which are simpler before proceeding to the more demanding terms 
(ii) and (iii). 

5.3.1 The term (i) 

The term (i) is equal to 



2 Jo 



E" 



a„ A"([/;,e:) / (z"(e) - z"(e:),$)^M"((c/r,e:),de) h{z:) 



ds 



^To be precise, here Dynkin's formula is applied to the map 

{U,e,u,p) ^ exp (i^^a^{U -u,ip)v +i^^{^,z"{0) -p)e] 

which is in the domain of the extended generator of the PDMP [U]^, 0", "(t), p(t))t>o satisfying the special condi- 
tions in Theorem 1 2. II 
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which can also be written as 



1 



E" 



/i(zr)(G"(c/r,e';)a>,<i>)£ 



ds. 



This term converges to 



where E[/i(Zs)] = (^(8,-0,$) is the characteristic function of the hmit process, due to 



- / {G{u{s),p{s))^,<i>)£E[h{Zs)]ds, 



a„E" 
t 



+ 



hiz:){G-{u:,e:)^,^)s 



ds 



{G{u{s),p{s))<5,<i>)£E[h{Zs)]ds 



W 



\h{z:)\ |a„(G"([/r,e:)<i>,$)£- (GKs),p(s))$,$)j 



ds 



{G{u{s),p{s))<^,^)£E[h{Z:)]ds- [ {G{u{s),p{s))^,^)£E[h{Z,)]ds 
Jo 



Indeed, the first integral in the right hand side converges to zero due to |/i(Z")| < 1 and condition 
(jC5.1|) . The second integral converges to zero as well due to the Bounded Convergence Theorem 
as E[/i(Z")] -^ E[h{Zs)] for all s > by definition of the weak convergence. 

5.3.2 The term (iv) 

The term (iv) is equal to 



^^/5^E" 



[A{z^{-), ^)s] {u:, e:) - (F(c/;, z"(e:)), $) J h{z] 



ds. 



It vanishes in the limit due to condition (jC5.2p and the following estimate: 



(iv) < 



^E" 
t 



[^(z"(.), $)£] (f/r, e:) - (F([/,r, z-i&:)), $), \hiz'^)\ 



ds 



< v^ E-\[A{z-{-),'^)s]{u:\e:)-{F{u:\z-{e:)),<^)^\ds. 







5.3.3 The term (v) 

The term (v) is equal to 






7^(z"(o - z"(e:), $)£) /."((c/r, er),do h{z:) 



In order to study its limit, note first that assumption (|C6.1[) implies the existence of a sequence 
/3„, n G N, converging to zerco such that 



lim E" 



A"(c/r,e^) 







'Sr|(<E',2:"(e)-2;"(ej)>£|>/J„ 



/i" (([/;, e:),de)ds 



= 0. 



(5.4) 



^For example, wc can construct such a sequence in the following way: Let a„(/J) denote the sequence elements in 
the limit of I ICS.lll which satisfy lim,i_>oo an(/3) = for all /3 > 0. The aim is to construct a sequence /3„ satisfying 
(A) lim„^oo ^n = and (B) lim„^oo a„(/3„) = 0. 

We first define a subsequence of /3„ by /3„^ = ra^^ for rra 6 N where the strictly increasing values of the indices 



Um has to be defined. For rim < n < rim+i wc set /3n 



This sequence converges to zero and hence 



the condition (A) on the sequence is satisfied. Choosing cleverly the values of the indices rim allows for also the 
condition (B) to be satisfied. 

We define these values recursively: Choose n\ to be the smallest integer n' such that a„(l) < 1 for all n > n' 
which exists due to IICS.lll . Then for the finitely many n < rii we can set /3„ arbitrary, say, equal to 1, as these 
sequence elements are insignificant for the limit. Now we have m = 1. Next set for rim+i some integer n" larger 
than rim which satisfies a„{(r?i + 1)~^) < {m + 1)~^ for all n > n" . For the n in between rim and rim+i we set 
0ri = m~^ and it thus holds that a„(/3„) < /J„ for all these n. The sequence a„(l3„) we construct in this way is 
monotonically decreasing and convergent to zero. 
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We then obtain for (v) the estimate 



E" 



t 



ds 



< fc / E" 

"'0 



A"([/:,ei') / aj(z"(0-^"(e:),$)£|V"((f/r,er) 

iv^|(2"(0-z"(ey),*)£|>/3„ 



ds 



+ ifc(/3„)|a„/ E" A"([/:,e^) / i(z"(o-z"(e:),$)£iV"((f/r,e:) 

Jo L "'v^f;ri(z"(?)-z"(e?),*>£i</3r. 

where k{u) := K{u)lv? is bounded by some constant k. The first term in the right hand side of 
this inequality converges to zero due to conditions (|C6.2[) and (j5.4p . as does the second term due 
to |fc(/3„)| = 0{l3n) in combination with condition (jC3.1|) . Therefore, the term (v) vanishes for 

n — > oo. 

5.3.4 The terms (ii) and (iii) 

The terms (ii) and (iii) are given by 



ds, 



u/aT.K'- 



{A{z"{Q:)) U: + BiUl\ z"{&:)) - Aip{s)) ^(s) ^ B(u(s),p(s)), i;)^ h{Z:) 



ds. 



which is further spht into parts containing the operators A and S, and 



i^c^E" 



(F([/;, z"(e:)) - F(«(s),p(s)), $), h{z:) 



ds. 



respectively. Let us start with the term (ii) containing and prove the convergence 



lim / iJ^W 

n-)-oo_/g 



(i?([/:,z"(e:))-B(u(s),p(s)),^)^M^r 



*(||,B:Ks),p(s)]^)^ + (||,i?;Ks),p(s)]V')^ds, 



(5.5) 



where the partial derivatives of ((5 are evaluated at {s^il^,^). Note that the operators i3*[u(s),p(s)] : 
X ^ X and B*[u{s),p{s)] : X ^ 8 are the adjoint operators to the extensions of the partial 
derivatives, see (CE]). 

The derivation of the limit for the operators A and F work completely analogously and the same 
line of argument yields the convergence to a limit given by analogous formulae as for the term 
concerning B, i.e.. 



lim / iV^E" (A(z"(eiO)C/r-A(p(s))7.(s),V)^/i(Z:') 



ds 



and 



(ff'^*(^(^)Wx + (|l'^^["(')'^(-^)]^>/-^ (^-^^ 



lim / ix/o^E'' 

n— >-cx3 



{F{Ul\z-{Q:))-F{u{s),p{s)),<^),h{ZT. 



As 



(5.7) 



{%KHs).v{sm) + 1%f;HsIp{s)] 



\d^j 



1$) ds. 
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where the partial derivatives of Lp are evaluated at (s, ?/>, $). Here the operators marked * are 
the adjoints to the extensions A{p{s)) e L{X*,X*), Ap[u{s) , p{s)] € L{£*,X*), Fu[u{s) , p{s)] G 
L{X*,£*) and Fp\u{s) , p{s)] <E L{£*,£*) of the partial derivatives. 

We now prove in detail (|5.5p . We first derive the pointwise limits for (almost) all s G [0,i] of the 
integrands in (|5.5p and then use the Dominated Convergence Theorem to infer the convergence of 
the integrals. Expanding the mapping i?(t/", z"(0")) around (w(s),p(s)) we obtain 

{BiU:,z"ie:))~B{u{s),pis)),^)^ =. (5.8) 

+ {B,Ms),p{s)]{U2~u{s)),i,)x + {B.p[u{s),p{s)]{z'\e':)-p{s)),i'}x+eB,n. 

We thus obtain - omitting the remainder term eB,n~ the equalities 

z^^,{B^[u{s),p{s)]{U:~u{s)),i:)_^h{Z:) ^ i^^,{B:[u{s),p{s)]{U^ ~ u{s)),^) ^ h{Z:) 

= tV^,{u:'-u{s),B:[u{s),p{s)]^)^h{Zl') 
/dh{Z^) 



{-^^B:[u{slp{s)]i)^eC 



and analogously 



zV^(i3,[w(s),p(s)](z"(e^)-p(.)),v^)^./i(z:) = (^MM,i?;[«(s),p(s)]^)^ e c. 

Here the operator _B*[u(s),p(s)] e L{X,X) is the adjoint of Bu[u{s)^p{s)] G L{X*,X*) and 
B*[u{s),p{s)] e L{X,£) is the adjoint oi Bp[u{s),p{s)] € L{£*,X*). Then the Continuous Map- 
ping Theorem and the continuity of the derivatives of h imply that 



(^,g:K.),pww)^ ^ (^,b:k.).pwi*). 



and 

where — ^ denotes convergence in distribution. Furthermore, the boundedness of h, the uniform 
integrability of Z" due to Lemma 14.11 and the continuity of the partial derivatives imply that the 
family of random variables given by the real and imaginary parts of the random variables 

{^^^,B:[uis),p{sM)^ , (^^!^,b;[u{s),p{sM)^ Ws<t,neN 

are uniformly intcgrable. Thus [HJ Appendix, Prop. 2.3] yields that 

Finally, the Dominated Convergence Theorem (see Appendix IA.3I for a detailed argument) allows 
to interchange taking the expectation and the Frechet derivative due to (j5.1[) and the boundedness 
of the exponential function along the imaginary axis. This yields 



e(^!^,b:[u{s)Msm)^ = {e^^,b:hs),p{sm)^ = (||,i?:[.(.),p(.)]^)^ 
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where the characteristic function of the hniiting process ip is evaluated at (s, ?/>, $). This holds for 
(almost) all s G [0,t]. Moreover, the integrands in the left hand side of ()5.5p are bounded in s 
and n, thus dominated convergence yields that the integrals converge and (|5.5p holds. We have 
deferred the detailed estimates of this last step to Appendix IA.3I 

Remark 5.2. Note that in parts (ii) and (iii) it was tacitly assumed that the remainder terms 
iA,n, iB,n and eF.ri hi the above expansions vanish for n — > oo. An argument why this holds under 
the Assumptions [C] is detailed in Appendix |B] 

5.3.5 Combination of the estimates 

We now study Lp'^{t,ip,^) — ip"{0,ip,^) which equals the sum of the previously discussed five 
terms. On the one hand, thanks to the convergence in distribution of Z" to a limit, the sequence 
(/3"(i, ip, $) converges to (p{t, ip, ^) for all i > and on the other hand the five terms converge to 
their limits respectively. Hence, we have obtained that the characteristic function (/? : M x X x £ ^ 
C : {t, ip, $) — >■ (p{t, ip, $) of the limit Z satisfies the complex partial differential equation 



= -i(G(7.(s),p(s))$, $)^ ^ + (^^,A*{p{s)) ^ + B:[u{s),pis)] ^ + F: Hs),p{s)] $) , 



+ (||, a;[u(s),p(s)] ^ + b;[u{s),p{s)] i, + f;[u{s),p{s)] $)^ (5.9) 

with initial condition (p(0, ■0, *&), where f and its partial derivatives -^, -g§- and -^ are evaluated 
at (i, i}}, $) and hence are elements of C, L(R, C) = C, L{X, C) and L{£, C), respectively. But this 
equation is precisely (j4.4p . Furthermore, as the equation (|5.9|) is linear in (p uniqueness of solutions 
to given initial conditions follows from the uniqueness of solutions started at (p{Q, ip, $) = 0. How- 
ever, the uniqueness of this solution is easily obtained by adapting the method of characteristics 
to partial differential equations in Hilbert spaces. We can thus uniquely characterise the one- 
dimensional distributions of the limit via the characteristic functions given by the unique solution 
to (|5.9p started at the limit of the initial distributions. At the end of this section, see Section [5. 4. 11 
we prove that (j5.9p uniquely characterises all finite-dimensional distributions. Let us consider first 
as an intermediate step a particular family of solutions to ((52 



5.4 Gaussian solution to the characteristic equation (|4.4p 

In this part of the proof we show that the limit possesses one-dimensional Gaussian marginals 
if the initial condition is Gaussian. Thus we recall that if a process Z is Gaussian, then its 
one-dimensional characteristic functions satisfy 

^(i, ^, $) = exp(i (mit), (iP, <i>))xx£ - i(C(i)(V', $), (^, $))xx£) (5.10) 

for families of m{-) E X* x £* and self-adjoint, non-negative trace class operators C{t) : X x f — >■ 
X* X £* . That is. (jS.lOp states that the one-dimensional distributions are Gaussian with mean 
m{-) and covariance C(-) [S]. We prove in the following that equation (|4.4p possesses a solution 
of the form (|5.10p . Note that we assume the initial condition is Gaussian, hence there exists a 
non-negative, self-adjoint trace class operator TZq and an element toq such that (jS.lOp holds for 
t = with C(0) = TZo and m(0) = mo. From here on we consider the spaces X, £ to be real 
Hilbert spaces. 

In order to obtain a solution of the form (jS.lOp for the equation (|4.4p (or (|5.9p ) we consider the 
ansatz 

(^(f,V,«') = exp{i(mi(t),V')x+i(m2(i),«')£ (5.11) 

-^[{Rlim,i')x + {Rl2m,^)x + (i?2l(t)VA$)£ + (i?22(t)$,$)£)} 
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where mi{t) £ X* and m2(t) G £* and the trace class operators i?ii(t) G Li{X,X*), i?22(i) € 
Li{£,£*), i?i2(t) e Li(£,X*) and i?2i(i) G Li{X,£*) are Frechet differentiable with respect 
to t. Note, that the adjoint of a trace class operator is again of trace class. Moreover, we 
can use these operators to define a trace class operator TZ E Li(^X x £,X* x £*^ such that 
S{t) :== (7^(t) (V'2, $2), W'l, *i))xx£ satisfies 

S{t) = (i?ll(i)^2,Vjl>X + (i?22W$2,$l)£ 

+ {Rl2{t)<^2,i'l)x + (i?2l(i)V'2, $l)£ (5.12) 

= {Rl,{t)i;i,ij2)x + {R*22mi.'^2)£ 

+ {RUt)iju^2)£ + {R*2imi,i'2)x ■ (5.13) 

The operator TZ is self-adjoint if the operators Ru and -R22 are and if R12 = ^21- However, 
these conditions are not assumed a-priori. Nevertheless, using (j5.12p we find that the second part 
in the exponent in (|5.11|) equals -'^(TZ{t) {tp, $), (?/', ^)) v- c ■ The strategy is now to show that 
there exists a solution to (|4.4p of the form (|5.1ip and that TZ{t) is a suitable Gaussian covariance 
operator C(i), i.e., 7^(t) is non-negative and self-adjoint. 

By the chain rule we obtain the partial derivatives of (p with respect to t,t,ip and $. Using the 
equality of the exponent in (jS.lip to (j5.13p we insert these partial derivatives into equation (j5.9p 
and compare the coefficients. We obtain for the mean values 

+ (m2(t), A;[u{s),p{s)] iP + B;[u{s),p{s)] iP)^ (5.14a) 

<^^^^)x = {m,{t),F:[u{s),p{s)]^)^ + {m2{t),F;[u{s),p{s)]^)^ (5.14b) 
and for the covariance operators 

{^^f^i^^i')x = {{Rii{t)+Rn{t))^,A*{p{s))i; + B:[u{s),p{s)]^)^ 

+ {{Rl2{t) + R2iit))^,A;[uis),pis)]^ + B;[u{s),p{s)]i^)^i5A4c) 

^dR^^^^^^ = (G(«(.),p(.))<i>,$)^ + ((i?i2(i) + i?;i(t))$,^:K5),p(.)]$)^ 

+ {{R22it)+R*22it))^,F;Hs),pis)]^)^ (5.14d) 

\ dt ' /^' 
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+ {R22m,A;[uis),pis)]^ + B;[uis),pis)]'4j)^ (5.14e) 

+ (i?*2(i)^,i?;[u(s),p(s)]$)^ 



(^*.*> 



X 



and 



+{R;2m, l;[w(s),p(s)]v + -§; Hs),p{s)] ^)^ (s.wf) 

+ (i?2i(t)V',i^;[^(s),p(s)]$)^ 

These real equations (|5.14p are satisfied for all tp, $ if and only if the mean values satisfy the 
following system of Hilbert- valued differential equations 

mi = {A + Bu)mi + {Ap + Bp)m2, 

(5.15) 
m2 = F„ mi + Fp TO2 

and the covariance operators satisfy the following system of operator- valued differential equations 

i?ii = i?ii o {A* + B*j + i?i2 o (I; + b;) + (A + Bu) o i?ii + {Ap + Bp) o i?2i, 

R22 = R22°F*+R2ioF*+FpoR22 + FuoRi2 + G, 

^ ^ ^ ^ ^ (5-16) 

R12 = Rii o F* + R12 o f; + {A + Bu) o i?i2 + {Ap + Bp) o i?22, 

R21 - R21 o {A* + B:) + R22 o (I; + B;) + Fu o Rn + Fp o i?2i . 

Note, that in the equations (|5.15p and (|5.16p we omitted the time-dependence of the derivatives 
in the right hand sides. Thus the systems (|5.15p and (|5.16p decouple in two closed systems of 
differential equations, the first takes values in the space X* x £* and the second is operator- 
valued. 

We next show that these systems possess a unique global solution to any given initial condition as 
non- autonomous (the derivatives being evaluated along (u(s),p(s)) depend on time) differential 
equations in Banach spaces. First, this is straight forward for system (|5.15p : by assumption the 
operators in its right hand side are linear and bounded. Hence, its right hand side is Lipschitz 
continuous. Moreover, the Lipschitz constant is bounded over [0,T] for every T >0, since the 
derivatives of A, B and F are continuous in t due to assumption (€[7]). Thus the Picard-Lindelof 
Theorem for Banach- valued differential equations, see, e.g., [11] Sec. 1.1], implies existence of a 
unique global solution of (jS.lSp to any initial condition. 

We proceed to system (|5.16p . The composition of a trace class operator with a hnear, bounded 
operator is again of trace class. Therefore we can interpret the composition of a linear mapping 
between spaces of trace class operators in the following sense. Let Yi, i = 1,2, 3, be Hilbert spaces, 
Oi € Li (1^1,12) be a trace class operator and O2 S ^(12,^3) a linear, bounded operator. Then 
Oi o O2 G Li(Yi, Y3) is a trace class operator and O20 : Li(Yi, ^2) — > -^1(^1, ^3) ■ O 1-^ O2 o O is 
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a linear map between two Banach spaces. Moreover, this map is bounded: 

< sup \\02\\LiY2.Ys)\\0\\L^iY^,Y2) = II 02!| L(yi,y2) 



Analogously, we obtain that for Oi G ^1(^27^3) and O2 G L{Yi,Y2), Oi 1— >■ Oi o O2 defines a 
linear, bounded mapping from ^1(12,^3) into Li{Yi,Y3). Therefore the set of operator-valued 
differential equations (|5.16|) is a linear, non- autonomous and inhomogcneous system (because of 
the term G{u{s),p(s))) in a vector Banach space of trace class operators. Thus it can also be 
analysed using the theory of ordinary differential equations in Banach spaces. The right hand side 
of the system (|5.15p . (j5.16p is linear, hence it is Lipschitz continuous with a uniform Lipschitz 
constant on any bounded interval [0,T] due to assumption {C^. As before the Picard-Lindelof 
Theorem yields existence of a unique global solution of (j5.16p to any initial condition. 

Finally, in order for the trace class operators Tl{t), t > 0, ~ the solution of the system (|5.16p with 
initial condition TZ{Q) = TZq - to define covariance operators of a Gaussian random variables we 
need additionally that they are (i) self-adjoint and (ii) that the map {ip, $) H' (7?.(t) (-0, <&), {ip, $)) 
is non-negative for all t > 0. By assumption the initial condition of the limiting process is Gaus- 
sian. Hence there exists a mean value mg G X* x £* and a non-negative, self-adjoint trace 
class operator TZq as its variance. Furthermore TZq is of the form (|5.12p for self-adjoint trace 
class operators i??^ e Li{X,X*) and i?22 G Li{8,£*) and operators i?i2(0) S Li{£,X*) e and 
-R2i(0) S Li{X^£*) such that R\2{^) = -R2i(0). The mean values and covariance operators of the 
one-dimensional marginals arc thus given by the solution to (|5.15D and (|5.16p . respectively, started 
at these initial conditions. Using the adjoints in the ansatz (|5.1ip we obtain the equations 

i?ti = Rii o {A* + Bl) + i?i2 o (I; + b;) + [A + Bu) o Rn + {Ap + Bp) o R21, 

R*22 = R220F;+R2l0F*+Fp0R22+Fu0Ri2+G, 



R*2i = Rn°F*+R,2oF; + iA + Bu)oR,2 + iAp + Bp)oR. 



(5.17) 



22, 



RI2 = R21 o {A* + Bl) + R22 o {A*p + B;) + F^ o i?n + Fp o R^i 



Comparing the right hand sides of these equations and of their respective adjoint operators in 
(j5.16p we find that they coincide. Thus we infer that Rii{t) and i?22(i) are self-adjoint and 
Ri2{t) ~ -^21(0 for any time t > as this holds for the initial conditions. Thus it remains to 
establish the non-negativity of the operator TZ{t). But the map exp{—^(TZ{t) {ip, $), {ifj, ^)) ^ A 
is the unique solution to (j5.9p with mean value function m{t) = and initial covariance TZq. 
Therefore it is necessarily the characteristic function of a random variable and consequently 
exp{— i(7^(i) {ip, $), {il>, ^)) xxE^ - ^' I^™n^ t'^is we infer that TZ{t) is positive. 

5.4.1 Conclusion of the proof 

Clearly, for a stochastic process to be Gaussian it is not sufficient that its one-dimensional distribu- 
tions, i.e., ()5.10p . are Gaussian but all finite-dimensional distributions have to be jointly Gaussian. 
Analogously, in order to uniquely define an arbitrary limiting process we have to characterise all 
its finite-dimensional distributions. Thus, we now complete the proof of Theorem 14 . 1 1 showing that 
all the finite-dimensional distributions are characterised via the solution of (j5.9p . We detail the 
step obtaining the two-dimensional from the one-dimensional characteristic function. The general 
finite-dimensional case follows by induction. 

Due to the Markov property of the PDMPs (t/", 0", u,p) or the processes {Z'^,u,p), respectively, 
the characteristic function of the two-dimensional distribution ip2 of Z" at fixed times ti < t2 
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satisfies 



(p^(tl,t2,V'l,V'2,$l,*2) 






:= E 

where E^ denotes the expectation with respect to the law of a PDMP Z" started almost surely in 
the point as indicated by the index. Clearly, this starting point has to be an attainable value for 
the random variable Z" . These are not necessarily all elements of X* x £* and we denote this set 
of values by /C„. We note that the function 5" thus defined is bounded and measurable. Next, we 
denote by g the characteristic function of -^ta-ti started almost surely at z and thus given by the 
solution to (|5.9|) at time ^2 ~ ii with respect to a degenerate Gaussian initial condition centred at 
z. The difficulty for taking the limit is that g" is usually not defined on the same set as g but only 



on a set K-n Q \/a{H x z"-{K"-)) possessing fuU measure, i.e. 
first extend g" to a function p" on all of X* x £* via 



[Z" e /C„] = 1. We thus have to 



9n{z) 



gn{z) ifze/C„, 
g{z) \i z ^ K-n 



With this extension it holds 



^^(tl,t2>l>2,$l,$2)=E" /l(Z,"^,Vl,$l)5n(^t","(il),P(il),i2-il,V'2,$2) 



Then we define the set 



5 := {z e X* X f * : g„(z„) -^ g{z) Vz„ -^ z} . 



Due to the foregoing convergence result and the continuity of g we immediately find that P[2^ti € 
5] = 1. Here we have essentially two cases for z" —> z. If all but finitely many elements of the 
sequence arc such that z" € /C„ then (7„(z„) — >■ g{z) holds due to the foregoing convergence result. 
On the other hand, if all but finitely many elements are such that z" ^ /C„ then g„(z") = g{z^) — >■ 
g{z) holds due to the continuity of g. All other possibilities are mixtures of such sequences and 
the convergence 5„(z") — >■ g{z) thus also holds. 

Therefore, [H Thm. 5.5] establishes that gniZ"^) converges weakly to g{Zt^) which yields for the 
two-dimensional characteristic function the convergence 

</?2(ii,i2,V'i,V'2,$i,$2) = fim «'[h{Zl,il;i,^i)g^{Zl,u{ti),p{ti)M - ti.^2,^2)\ 

= E[/i(Zt,,Vi,$i).9(^ti,u(ti),p(ti),i2-ii,V'2,$2)] . 

Here the left hand side Lp2 is the two-dimensional characteristic function of the limiting process Z . 
We conclude the proof showing that v?2 is uniquely defined from (|5.9p . We have that g{z,ti,t2 — 
ti) ~ omitting the remaining arguments - is the characteristic function given by the solution of 
(|5.9p started at the characteristic function of the degenerate Gaussian law centered at z, i.e., 
the initial condition to (j5.9|) is given by exp(i(('02, ^2), -z)), at time ti in z after time ^2 ~ ii- 
Thus g{z,ti,t2 — ti) is of the form (|5.10p where the mean value function at time ^2 ~ ti is given 
by the solution to (|5.15p started in ti at z. Due to the general existence theorem of linear 
equations in Banach spaces there exists a linear, bounded operator M{ti,t2) such that this mean 
function is given by M{ti,t2)z, i.e., the first part of the exponent is (M(ti, t2)-^ti, {i^2,^2))xx£ = 
{Zt-^, M{ti,t2)*{ip2,^2))xx£- Secondly, the covariance operator 7?. at time ^2 — ii is given by 
the solution to (|5.16p started in the origin at ti which we denote by TZo{ti,t2)- Note that it is 
independent of z. Therefore we obtain 



'/52(il,i2,V'l,V'2,$l,$2) = E 



^i(Ztj,(l^l,*l)4-A/'(tl,t2)(l/'2,*2))xx£ 

X exp 



(5.18) 

{-i(7^o(tl,i2)(^2,$2),W'2,$2))xx£}, 
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where the remaining expectation in the right hand side is a one-dimensional characteristic function 
uniquely given by (|5.9|) . Hence, all two-dimensional characteristic functions are uniquely defined. 
The analogous result for all finite-dimensional equations follows by induction. We furthermore infer 
from ()5.18|) that the limiting process is a Gaussian process, if the limit of the initial conditions is 
Gaussian. 

5.5 Proof of Corollary 14.11 

If the initial conditions of the PDMPs admit a Gaussian limit then the limiting process is under 
Theorem 14.11 a Gaussian process as is the mild solution to an stochastic evolution equation of 
the form (|4.6|) started in a Gaussian initial condition (cf. [3]). Hence, it suffices to show that 
the finite-dimensional characteristic functions coincide. In the following we first show that the 
one-dimensional marginals coincide, i.e., we show that the mean and covariance of the solution of 
the SPDE satisfy (j5.15|) and (j5.16|) . respectively, and afterwards extend the result by induction 
to all finite-dimensional distributions. For future use we repeat these equations which we write in 
simplified matrix notation, i.e., for the equation of the mean values (jS.lSp we write 

mi \ /I+B„ Ip + BpWmiX 

where the matrix-vector multiplication is understood as applications of operators, and similarly 
for the equation for the covariance operator (|5.16l) 



i?ll 


R12 


] 






R21 


R22 


1 






= 


fRii 

\R21 


R12 \ 
R22 J 


\a; + b; 


K 

F* 
p 



(5.20) 




Bp \ ( Rii i?i2 V /O 



Fp / v^2i i?22 ; VO G 



where the matrix-matrix multiplication is understood as composition of operators. We denote by 
S{s, t), s < t, the two-parameter semigroup of linear operators on X* x £* defined by the solution 
of the differential equation for the mean, i.e., S{s,t)m is the solution to (j5.19p at time t started 
in m in time s. Then the unique mild solution to (j4.6|) is given by 




S{s,t)l ]dW^, (5.21) 

where W'^ is a Wiener process on a Hilbert space H' with covariance Q G L{H'), g{t) E L{H' , £*) 
and denotes the origin in L{H' ,X*). Thus taking expectations on both sides of (|5.2ip and 
differentiating with respect to t we immediately obtain that the mean value function of the mild 
solution (|5.2ip satisfies (|5.19p with initial condition given by the expectation of (f/p, Pq). 

It remains to consider the covariance operator. We denote by Qo the covariance operator of the 
initial condition and by Qy (t) the covariance of the stochastic convolution 

Yt^ f Sis,t){0,gis)fdW? 
Jo 

which is given by (cf. [5]) 

QY{t)= I (5(i,s)(0,5(s))^)g(5(i,s)(0,5(s))^)*ds eL{X*x£*,X*x£*). 

Thus the covariance of Zt is given by Qz{t) = 5(0, i)(5o'5'*(0, i) + Qvit)- Here all the adjoints are 
understood as Hilbert space adjoints. Then Qz{t) satisfies the initial condition Qz(0) = Qoi and 
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differentiating with respect to t yields 
dQzjt) _ ( A + Bu Ap + Bp 




Qz{t)+Qz{t)r7.^~: ~: + n n ■ (5.22) 



dt \ F„, Fp 

Thus Qz satisfies the differential equation for the covariancc of the limit (j5.20p . Thus the mild 
solution of the SPDE possesses the same one-dimensional marginals as the limiting process char- 
acterised by (j4.4p . The Markov property of the SPDE solution yields that the two-dimensional 
characteristic functions are equal to (jS.lSp . The general statement that the finite-dimensional 
characteristic functions agree with the finite-dimensional characteristic functions of the limit char- 
acterised by (|4.4p follows by induction. The proof is complete. 

6 Applications to Electrophysiology 

In this section we apply the main result, Theorem 14.11 to obtain limits for the global fluctua- 
tions of stochastic microscopic models around their deterministic limit. These limits are classical 
deterministic spatio-temporal models in neuroscience. The first one is a simple version of the 
Hodgkin-Huxley model (cf. [TS]) for an excitable membrane reduced to its essentials, which is 
used not only to model neuronal membranes but also intracellular tissues in calcium dynam- 
ics or cardiac tissues. Secondly, we consider a neural field model for the large scale activity in 
macroscopic brain areas. The corresponding deterministic limit in this case is the Wilson- Cowan 
equation [7] . Internal fiuctuations for such models have been considered in previous work applying 
the law of large numbers (cf. Theorem 13. ip and the martingale limit theorem [571 US] • For these 
models the spatial domain is a bounded domain D CM.'^ with sufficiently regular boundary. The 
Hilbert spaces that thus occur as state spaces of the PDMPs are spaces of real functions on D, 
for instance, the Lebesgue space L'^iD) or Sobolev Spaces H°'{D). Subsequently, with regards to 
brevity and simplicity of notation we simply write L^ or 7?" for these spaces, omitting the domain. 
Furthermore, products of elements of these spaces are always understood as pointwise products of 
real functions. 

Remark 6.1 (Spatial discretisation). In both applications we will study, we start with discretising 
the spatial domain D. For each n <S N we denote by V" a decomposition of D into p{n) < oo 
non-overlapping subdomains £)'''", k = l,...,p(n). For technical reasons we assume that the 
subdomains £)'''" are convex [5B1 [551 [57] . Then, i'±{n) and 6±{n) denote the maximum / minimum 
Lebesgue measure and the maximum / minimum diameter of the subdomains in partition P", 
respectively. We refer to our previous work for a discussion of the existence of such partitions 
satisfying all the subsequent assumptions for a large class of domains containing all practically 
relevant geometries. 

Remark 6.2 (Hilbert scales). Another important concept are Hilbert scales. A Hilbert scale is 
a family of Hilbert spaces Ha indexed by a e R^ which satisfy Ha-^ ^-> iJ^a for cti > Q!2- This 
extends to Hilbert scale indexed on a G M if we set H^a the dual space of Ha with the following 
identification of inner product and duality pairing 

{(f>,x)H_^^ [4i,x)ho y (jye Ho, X e Ha- (6.1) 

As we observe from this identification the Hilbert space Hq plays a distinguished role as each 
other space Ha^ a > 0, forms an evolution triplet with i/o- Examples of Hilbert scales are the 
standard Sobolev spaces H°' , a € N, on a bounded domain D or Sobolev spaces Hq of functions 
vanishing at the boundary with the usual interpolation spaces for non-integer indices [T] . Here the 
distinguished Hilbert space is given by H^ = L^. A second example arc Hilbert scales obtained 
by closing of operators, see, e.g., [S] [121 HO] in applications to limit theorems for spatio-temporal 
reaction-diffusion models. 
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As an example consider the one-dimensional negative Laplacian —A with Dirichlet boundary 
conditions D = [0,1], and the spaces defined by 

oo 

H^ := {geL': ^^(g, ^,)h (1 + *')"} , (6.2) 

where ipi{x) = \/2/l sm{'!Tix / 1) , i > 1, form a complete orthonormal basis in L^ consisting of 
eigenvectors to the Laplacian. The inner product on this space is given by 

oo 
i=l 

Here, as before, the distinguished space satisfies Hq = L^. As the Sobolev spaces above, the spaces 
(|6.2p create a Hilbert scale. Moreover, these spaces are subspaces of the respective Sobolev spaces 
iJ", where the norm induced by (•, Off-c. i^ equivalent to the restriction of the usual Sobolev norm 
PP] , Thus in terms of embeddings and forming of evolution triplets this scale behaves as the 
standard Sobolev spaces. Another important result is the following [301 Lemma 2.2]: the Laplace 
operator A : H^ — > H~^ possesses an extension A_q, to an operator acting on H^a for every 
a> 0. 

Finally, we note that in the following we always use the notation Ha to denote the spaces defined 



6.1 Compartmental-type neuron models 

Models of excitable membranes describe the evolution of the transmembrane potential, that is the 
difference 14xt — ^nt of electrical potential between the outside and the inside of the cell. This 
evolution results from a complex electrochemical coupling between the membrane potential and 
the various ion channels located at some points of the membrane. The repartition of these channels 
is heterogeneous. In particular for neuronal membranes it is known that channels are present only 
at the Nodes of Ranvier. We first present a simple version of general excitable media models. 
However, the results extend immediately to more general versions (cf. [25l[27])- For the spatial 
domain we restrict ourselves to a finite interval D = [0,1] C M. For this model the deterministic 
limit is given by the excitable membrane equation 

it ~ Au + p ■ (v — u) 

(6.3) 
p = a{u) • (1 — p) — b{u) ■ p 

on [0,1], with Dirichlet boundary conditions to conform with the exposition in Remark 16.21 In 
(|6.3p A denotes the Laplace operator, TJ is some constant and a,b : M. —>■ M_|_ are non-negative 
smooth functions. The model (|6.3p is stripped of any complicating constants and corresponds to 
an excitable membrane with currents due to a single family of ion channels which can be in two 
states, i.e., they are either open or closed. It was shown that systems of the form (j6.3p possess 
unique global solutions which are componentwise in C{W+,H2) and C{M.^,H'^), respectively, for 
all initial condition uq G H2 and po G H^ f5S] . 

In the corresponding stochastic PDMP model, we get the family of abstract evolution equations 

ii = Au + 2"(r).(w-u) (6.4) 



equipped with the same boundary conditions as the deterministic model (|6.3p . This yields that 
we choose for X d H d X* a standard evolution triplet in order to deal with reaction-diffusion 
equations, i.e.. 
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Note that Maurin's Theorem imphes that the embedding Hi ^^ L^ is of Hilbert-Schmidt type 
as the spatial dimension is one. The piecewise constant component d" ~ (0^'", . . . , ^p*^")'") S 
Kn C NP(") counts the number of open ion channels in the subdomains Z?'^'-". The Central Limit 
Theorem that we proved in the previous section enables us to specify the size of the intrinsic noise 
induced by the stochastic gating of the channels using parameters of the domain decomposition. 
One interest of our result is that we can, up to some extent, consider inhomogeneous channel 
repartition. More details are given in Theorem 16. II below. 



Given a partition P" of subdomains D*-''", let us denote by l{k,n) < od the total number of 
channels in domain Z?*^'". We denote by i±{n) the maximum / minimum number of channels in 
the individual subdomains of V". The transition dynamics of the piecewise constant components 
0" are given by the (M,0")-dependent rates 



0"- with rate < 



a(u'^^") {l{k,n)- 




if ^^" = e''-" - 1, e^-" 


= 03-", J ^ k, 


if ^^" = e'''" + 1, 03'" 


= 0^'", J ^ k, 


otherwise , 





where u '" = |D''''"| ^ /n&.n u{x) dx is the average membrane potential over the fc-th subdomain 
of P". We see that the components of 6" satisfy 6''^''" G {0, . . . ,l{k,n)} and hence X" is finite. 
Finally, we repeat that the coordinate functions z" are given by (j2.5p . i.e., 

P(") nk.n 

and thus E ^ L^. It remains to fix the Hilbert spaces £ and V: we chose V = H" and £ = H^" 
for a e (1/2, 1] and Maurin's Theorem implies that the embedding £ ^^ y is of Hilbert-Schmidt 
type. 

If the initial conditions converge in probability and i{n)- — > oo, (S(?t.)+ — > the sequence of 
stochastic models (?7", 2;"'(9")) converges u.c.p. to the solution (u,p) of (|6.3p in L^ x L^ [57]. We 
now present the central limit theorem for the global fluctuations. 

Theorem 6.1. Let us assume that the rescaled difference in the macroscopic initial conditions 
converges to a Gaussian distribution Ao and that i-{n) — >■ oo, <5+(n) — > with £_(n)(5+(n) — > 

and lim„_>.oo f(n)^^(n) ~ ^- ^^^ '^ ^ (l/2j !]• Then the sequence 



converges weakly in iJ_i x H~^" to the mild solution of the SPDE system 



AUt = 



At 



(6.5a) 



{^^p{t))Ut + {v-u{s))Pt 

dPt = F'[uit),p{t)]iUuPt)dt +git)dWt^ 

with 

F'[u,p]{U,P) = {a'{u){l-p)-b'{u)p)U -{a{u) + b{u))P (6.5b) 

and initial conditions {Uo,Po) distributed according to JVq on iZ_i x H~^". The process W^ is 
a cylindrical Wiener process on a Hilbert space H' with covariance operator Q S L[H') such that 
g{t) E L{H' , H^^") satisfies g{t)Qg*{t) = G{u{t),p{t)) o L^^a, where G is defined via the bilinear 
form 

{Giu,p)i;,<l))H-2^ = na{u{x))il-p{x))+biu{x))pix)')ijix)(t)ix)dx. (6.5c) 



D 
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Note that the assumption that hm„^.oo ^ ["i^ M = 1 in the above theorem ahows for hetero- 
geneity in the spatial repartition of channels in the membrane. 

Proof. We have to establish conditions (CJT]) - {(^ for the choice of Hilbert spaces H = E = L^, 
X = Hi, £ = H^" and V = H". Then the convergence will follow from Theorem 14.11 and Corollary 
EH 

Some of these conditions are easier to check than others. Note that (C[2l) is satisfied by assumption. 
Moreover as the embeddings Hi ^-> L^ and H^" ^-> H" are of Hilbert-Schmidt type, condition 
(CSl) follows from {CQ - {C^ as discussed in Section O Conditions (|C5.1|) and (CE]) were 
established in [57] to which we refer the reader. The PDMPs as well as the deterministic limit are 
pointwise bounded over D independently of n which ensures condition (C|5|) holds. 
Hence, we are left with (CHJ, (CjS]), (|C5.2p and (CjT]) which we now establish. 

(a) Lipschitz conditions in (CHJ : From the coercivity of the Laplace operator we deduce 

{A{U'^-u),U--u)h^ < -Ji\\U--u\\l,+j2\\U--u\\l.. 
for some constants 71,72 > 0. Next, we obtain for 1/2 < a < 1 that 

(z"(e") {v - C/") -p{v-u),U'^^u)h^ - - (z"(e"), ([/" - uf)^., + (z"(e") -p,(v^ u){U- - u))h^ 

< \\z'^&")-p\\h-^\\v-u\\h^\\U^'^u\\h^ 

where C is a constant resulting from the continuous embedding H^ ^-> H". Note here, that for 
these estimates to hold we have to choose a < I. This restriction combined with the further 
condition a > d/2 where d is the dimension of the domain D, impose that we take d = 1. 
Integrating the above two estimates over [0, T] we find that if 

5 := -71 / ||C/"-u||^idi + eC / \\v - uW^r \\U" - uW^r dt < 
Jo Jo 

then the first Lipschitz condition in (Cnj is satisfied. Standard estimates in the theory of linear 
partial differential equations, see [13| . yield the bound 

l|w|U-((o.T),ffi) < C{\\uo\\%i +\\p- {v-u)\\l2((o,t),l'^)) < C{\\uo\\jji +Tv) . 
Hence \\v — u\\l°°{{o,t),h^} is finite and we estimate S by 

S<-li[ \\U"-u\\l^dt + eC\\v-u\\l^..,J..H^. f |it/"-«||2,,dt 
Jo Jo 

Now, choosing e sufficiently small we find that S < and thus there exists a constant L3 such 
that 

i-T 

(A([/" ~u) + z"(e") {v - [/") - p • (iJ - u), J7" - u)h^ dt 

< L3 ^||C/"-^.l|i. + |i^"(e")-p|l?,_„di. 
Jo 

Hence the first Lipschitz condition in (CH]) is satisfied. In order to derive the second one we use 
that ||z"(e")||L~ < 1 and obtain 

\\F{u,p)-F{U'\z^)\\h-^ 

< Wiaiu) + b{u))p- (aiW^) + biUn) ^"(e")lk-« + \Hu) - a(t/")||H- 

< II (a(«) + b{u)) {p - 2"(e"))||H-» + C II ia{u) - a{Un + Ku) - h{U'')) z^iQ^U^- + C \\a{u) - a(C/")|U2 

< ||(a(w) + b{u)) ■ {p - z"(e"))|U-» + iCKb \\u - [/"lU., 
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where C is a constant resulting from the continuous embedding of L^ into H " and Lab is a 
common Lipschitz constant for a and b. Finally, since 

||(a(ii) + 6H)(p-z"(e"))||H-. = sup \(v,{a{u) + h{u)){p~z-))^,\ 

||'Ii||h-q=1 

< sup (\\v{a{u) + h{u))\\H^\\p-z^H-A 

Ht'll/i-a— i 

< C\\a{u) + b{u)\\H^\\p-z^{Q'')\\H-'^, 

where C is a constant resulting from the continuous embedding H^ ^^ if", we find the second 
Lipschitz condition in {C^ is satisfied since ||a(u) + 6(M)||ioc((o,T),ffi) < oo- We estabhsh the third 
Lipschitz condition in (CJT]) by the same estimates with a replaced by 2a. 

(b) We now estabhsh (C^ and (fC5?2|) using that (see [37]) for a > 1/2 



E" 
/o 



A"(f/;,e:)/ |jz"(0-^"(e;^)||l,-„/i"((f/r,er),de)ds =o(^+(n)/£_(n)), (6.6) 



T 



and 

/ E"|| [^"(z;(-), • )l2] ((7", 9") - F(t/", z"(e"))||^, ds = 0(<5^(n)) . (6.7) 

JQ 

The asymptotic behaviour ()6.6[) implies (jC3.1|) . The asymptotic behaviour (j6.7p combined with 
the continuous embedding L^ M' iJ^" and the assumption ^_(n)(5+(ri) -^ imply (jC3.2p . The 
same reasoning implies condition (|C5.2p . 

(c) It remains the consider the differentiability conditions (CJT]). Remember that in this section 
the operators A, B and F are given by 

A{p) ~ A, B{u,p) — p{v — u), F{u,p) = —(a(u) + b{u)) p + a{u) . 

We have already stated that Laplacian can be extended to an operator A_a g L{H-a,H-a) 
(cf. Remark 16. 2p . As it is independent oi p, A = A_q, trivially satisfies the differentiability 
condition (Cj?!). It is easy to see that the operator B is partially differentiable with derivatives 
Bu[u,p] G L{Hi,H^i) and Bp[u,p] € L{L'',H^i): 

h I— > Bu[u,p\h = —ph, h I— > Bp[u,p\h = ( v — u)h . 

As p and v ^ u evaluated along the deterministic solution are elements of H^ at each time t, these 
operators immediately extend to operators Bu[u,p] S L(iJ_i,if_i) and Bp[u,p] £ L{H^^°',H^i) 
being densely defined and bounded. For instance, the derivative Bu[u,p\(j) — —p4> G ^-i for 
(j) e H-i is defined as 

{-p(f),v)Hi ■= {4),-pv)Hi \fv e Hi. 

This is well-defined as the pointwise product pv is an element of Hi due to the Banach algebra 
property of this space on one-dimensional domains. Like the operator B, the operator F is linear 
in p and the derivative Fp[u,p] e L{L^,L^) is given by /i i-> Fp[u,p]h — — (a(u) + b{u)) h which 
extends in the same way as Bp above to an operator Fp[u,p] E L{H^^",H^^°'). Dependence on 
u is nonlinear, thus the derivative Fu[u,p] £ L{H^,LF') is 

h^ Fu[u,p\h = -{a'{u) + b'{u))ph + a'{u) h 

which completely is extended to an operator Fu[u,p] in L{H^i,H^^") analogously. D 
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6.2 Neural field models 

As ill the previous eompartmeiital model, one chooses a sequence of partitions P" of a spatial 
domain D C M''. The sequence of PDMPs is given by a sequence of jump processes connected 
to this sequence of partitions. For this model there is no restriction on the dimension of the 
domain. The main difference currently is that here we only consider jumping components 0" 
and no continuous component [/". The deterministic limit is here given by the solution to the 
Wilson-Cowan model which describes the mean level of activity of populations of neurons and 
reads as follows: 

p{t,x) = -p(t,x) + fH w{x,y)p(t,y)dyj , (6.8) 

The gain function / and the connectivity kernel w are assumed to be smooth and bounded. As 
previously mentioned only finite number of neurons are accessible so it is natural to consider 
stochastic models in which the piecewise constant components 9" — (0^'", . . . , fl^'"''") G A'„ C 
pjp(n) count the number of active neurons in the individual subdomains. Here l{k, n) is a parameter 
which is related to the number of neurons in the subdomains in the sense that l(k, n) > E"0q'". 
The transition dynamics of the components 0" is given by 

p{n) 



6*" with rate < 



Y^ tU" \ ;f Tik.n _ nk.: 






otherwise , 



where 



Finally, as in the previous applications the coordinate functions are given by 

P(") Qk.n 

^"(0 = E^(^I..^.i^^ (6.9) 

The Hilbert spaces are chosen as before, i.e., E = L^, V = H" and £ = H^" with a > d/2. If the 
initial conditions converge in probability and £(n)^ — > cxj, 6(n)^ — >■ 0, the sequence of stochastic 
models z"(6") converges u.c.p. to the solution p of the the Wilson-Cowan equation (j6.8p in L^. 
This was proved in [55] building on [37] ■ We now present the central limit theorem for the global 
fluctuations. 

Theorem 6.2. Let a > d/2 and assume that the rescaled difference in the macroscopic ini- 
tial conditions converges to a Gaussian distribution Ao and that £-{n) — >■ oo, S^{n) — > with 
£-{n)6+{n) — >■ and lim„_j.oo ^~;'^^^^;"^ — 1. Then the sequence 



v1tS(^"(©")-^) 



converges weakly to an H '^"'-valued diffusion process which is a version of the mild solution to 
the SPDE 

dYt = Fp[p{t)]Yt At + g{t) AWu (6.10a) 

where Fp[p{t)] e L(iJ-2«^//-2a) jg defined by 



{Fp\p\h,v)H2^ = -(/i,w)h2» + / i^f'{{w{x,-),p)L2) {h,w{x,-))H2ojv{x)dx V/i,w 



-2q 



„ ^ .,„ ^ , , , . eH- 

(6.10b) 
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and initial conditions Pq are distributed according to Ao on H ^" . The process W^ is a cylindrical 
Wiener process on a Hilbert space H' with covariance operator Q G L{H') such that g{t) G 
L{H',H-^°') satisfies g{t)Qg*{t) = G{p{t)) o t'L, with G{p) € Li{H^" , H'"^") defined by ' 

(G{p)h,v)H2. = J (p{x) + f(^J w{x,y)piy)dy)^h{x)vix)dx yh,veH^^ (6.10c) 

Proof. As for the compartmental models we only need to establish (Cjl]), (03]), (jC5.2[) and (CJT]). 
For the other conditions see [5B]. We start with the Lipschitz conditions in (CH]). It is shown in 
pS] that the function F{p) given by the right hand side of (j6.8[) satisfies a Lipschitz condition in 
every space H~°', a > 0, i.e., there exist constants L_q such that 

\\F{gi,t) - F{g2,t)\\H-^ < £-a!|.gi - .galk-^ Vt > 0,51,^2 e L\ 

Therefore condition (C[T]) holds. The same asymptotic behaviours as (|6.6p and (|6.7p for the neural 
field case are proven in [55]. Analogous arguments yield again conditions {C^ and (|C5.2[) are 
satisfied. We now consider the differentiability of the right hand side F of (|6.8p to establish 
condition (CE]). It is straightforward to see that F is strongly Frechet-differentiable in L^ for 
/ G C^O^) with derivative Fp[p] € L{L^,L'^) given by 

h ^ Fp[p] h := -h + f n w{-,y)p{y) dyj / w{-,y)h{y)dy . 

For sufficiently smooth y h- !> w{x,y) this map extends to an operator in L{H^^°', H~^") via the 
definition 

{Fp[p]h,v)H2^ = -{h,v)H2.+ J [f'{iw{x,-),p)L2) {h,wix,-))H-^) v{x)dx . (6.11) 

Therefore condition (C[7]) is satisfied. D 

A Appendix 

In this Appendix we have gathered detailed results and remarks mentioned without details in the 
main text. 

A.l A function satisfying the assumptions of Theorem 12.11 

We show here that the function / given by (j2.4p satisfies the assumptions of Theorem 12.11 (b) . 
This function is complex which however, does not pose any difficulties as C can be identified with 
M^ and / considered componentwise in which case it possesses a bounded real and imaginary 
part as |/(u,6')|c ~ 1 for all {u,9) E H x K. As / is bounded the integrability with respect 
to the compensating measure is immediate. Moreover, its integral with respect to the associated 
fundamental martingale measure is a martingale due to [171 Thm. 4.6.1] and thus Dynkin's formula 
12.31 holds. Furthermore, / is partially continuously Frechet differentiable with respect to the first 
argument due to the chain rule and the Riesz Representation /„ of the partial derivative is given 

by 

fu{u, 0) = i exp(^i (^', u)h + i ("I", z{0))e) i^ e H* , 

where H is considered a complex Hilbert space. As ■;/) G A" it holds that fu{u,6) G X and finally 
it is easy to see that /„ is a locally bounded composition operator as 



iexp(i(V',w)H + i($,z(0))B)V < Pyjfx 

for aU {u,e) G i2((o,T),A X A^. Hence, formula 1^ holds for / as defined in ^^ 
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A. 2 An auxiliary convergence result 

Let us prove a convergence result for the trace of the quadratic variations of the PDMPs. This 
result was previously stated without proof in [27] . 

Assumptions B. (Bl) For all T > 0, 



sup a„ E" f Tr G"{U^, 9^) ds < oo , 

nGN Jo 



where G"'{u,9"') is the operator from £ to £* defined in (j2.7|) . 
(B2) There exists an orthonormal basis {(pk)kei'i of £ such that for all A: e N and all T > 



a„ E'' 



(G"(t/r,e:)(^fe,V:>fe)£ds < jkC{T). 



The constants 7fc > are independent of n, T and satisfy X^fceN^'^ ^ °*^- "^^^ constant 
C{T) > is independent of n, k. 

Note that assumption (mj is implied by (|C3.ip and assumption (E(5|) is in applications usually 
an intermediate result in establishing tightness of the PDMP sequence. We can then prove the 
following result. 

Proposition A.l. The Assumptions [Bl and the convergence (jC5.1[) imply 

lim a„ / E"TrG"([/;,e;')ds = / TrG{u{s),p{s))ds . 
"-*°° Jo Jo 

Proof. Let ($fc)fegN be an orthonormal basis of £ then it holds by Jensen's inequality that 



l(G(^i(s),p(s))*fc,$fc)£-«„(G"(c/r,e:)$fc,$fc)^|ds 



> 



(GKs),p(.s))$fc,$fc) d.s. 



'(G"(c/:\e;')$fc,$fc) ds 



As the left hand side converges to zero due to (|C5.1I) it follows that for all $fc, fc e N, 



lim a„ 

n— foo 



E"(G"(C/:,e:)$fc,<i>fc)^ds = / {G{uis),pis))^k,<^k),ds. 
/o Jo 

This implies that any finite sum of terms in the right hand side converges to the corresponding 
finite sum of the limits in the left hand side. Moreover due to dominated convergence, 

rT i-T 



j W'TvG'^{U^,Q1)As = / E"^(G"(t/;\e^)$fc,$fe)^ds 

= 5] /'E"(G"([/;,e:)$fc,$fc)^ds 



(A.1) 



and analogously 



/ TrG(u(.s),p(s))ds-;^ / (G(u(.s),p(s))$fe,$fc)^ds. (A.2) 

Jo fcgi^Jo 

We show now that the partial sums for the sequence of PDMP models converge to the trace 
uniformly in n: 

Ve>0 3A/ : Vm > Af : 



E 



an I E"(G"(t/:\e;')$fc,$fc) d.s-a„ / E"Tr G"(C/r, e^') ds 











< e Vn e N. 
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Note that (jA.lj) implies 



E 



k<r. 



E"(G"(c/;, e;^) $fe, $fe) ds - «„ / E"Tr g"(c/;, e;*) ds 



E- 



E"(G"(C/r,e^)$fe,$fc),ds. 



Due to condition (E(2|) there exists a sequence of 7^ > independent of t,n with J2keN^k < 00 
and a constant C{T) > independent of n such that 

^a,,! E"(G"(C/r,e:)$,,$,)^ds < G(T)^7fc. 

k>m fc>?n 

As X^feGN 7^ "^ ^^ f^^ every e > there exists an Af such that 

C(T)^7fe < e ym>M, 

k>rn 

where M does not depend on n e N. This yields the uniform convergence. 

Uniform convergence is important because it allows to change the order of taking limits. Thus we 
obtain 

hm hm V«„/ E"(G"(t/r,e:)$fc,$fc),ds 

k<m ''" 

= lim lim V«„/ E"(G"([/r,e:)$fc,$fc)^ds 

k<m -^0 

= hm ^ / (G(tx(s),p(s))<i>fc,<i>fc)^ds 

rn— >-cx3 ■^^ — ^ /n 



TrG(M(s),p(.s))ds, 



where we have used (|A.2[) for the last equality. The proof is completed. 

A. 3 Bounds for the application of the Dominated Convergence Theorem 



D 



We now turn to some detailed estimation procedures that we did not include in the proof of the 
Central Limit Theorem: bounds for the application of the Dominated Convergence Theorem, for, 
on the one hand, interchanging differentiation and taking expectation and, on the other hand, for 
concluding from pointwise convergence to convergence of the integrals. In the following H denotes 
appropriately either the space X or £. Analogously (j> denotes either ip £ X or ^ ^ £, respectively, 
and so does Z^ denote either the X*- or the £'*-component of the limiting process Z. 

We start presenting an argument which allows to interchange Frechet differentiation and taking 
the expectation. Assume that 



lim 

Ikllw-yO 



which is equivalent to 



lim 



^hjZs, <p + x)- h{Z,,4>) - {^^[<f\,x) 



ip{s, (t) + x)- ip{, s(j)) - (E^^^ [4>],x)^ 



(A.3) 



m\H 
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Thus the uniqueness of the derivative imphes that gri^] — E — nr^^i^] which is what we aim for. 
Hence, in order to obtain ()A.3p it suffices to show that 

\\h{Z,,(t) + x)-h{Zs,(j})-{^H^[4>],x) \\ 
lim Hi — = . 

Ibllw^n \\x\\h 

However, here the term inside the expectation converges to zero almost surely, hence we can 
employ the Dominated Convergence Theorem to conclude from pointwise limits to limits of the 
expectation. That is, we are left to find an integrable upper bound to the terms 



\h{Zs,cl) + x) - h{Zs,'iJj)\\c 



and 



dhiZs) , 



L{H,C) 



In order to bound the first term we employ the Mean Value Theorem and obtain 



\h{Z,,c^ + x)-h{Z, 



\h\ 



n 



< 



X\\-H 
1 



iZl'-"h{Zs,(j) + rx)dr\x 



\iZ^ h{Zs,ij + rx)\\ dr 



< \\z\ 



n'l 



For the second term we obtain 
dhiZs) 



[</>] 



«* 



Wi 



= \\ih{z,,c^)z^ U^<\\z^ U 



L(H,C) 

Clearly, the upper bound is integrable, thus we have completed the argument. 

For the second application of the Dominated Convergence Theorem in order to conclude from 
pointwise convergence of limn^oo E"^ „' ' ,g{s)'^ ~ E<^ — ai 'ff('^))-f/ ^'-"" almost all s < t to 



d^p 



dip 



we need the integrands in the left hand side to be bounded in n almost everywhere by a function 
integrable over (0,i). Here g{s) represents the points where the derivatives are evaluated such 
as, e.g., B'^[u{s),p{s)\(j). As these derivatives are continuous - as is the deterministic solution 
{u{s),p{s)) - it holds that sup^jgrp ^i ||.g(s)||-H < oo. Simple estimates yield 



^"(^^"^(^))« < ^r. 



dh{z^,^). 



dtp 



HA\a{s)\\n 



< sup ||.g(s)||«E"||*Zr'^ h{Z,,ij)\\H' 

sG[0,t] 

< sup |J3(s)l|«E"l|Z:^-«*|j«.. 
se[o,t] 

Here the last expectation in the right hand side possesses a uniform bound independent of n G N 
and s <t due to Lcmma l4.1l Therefore there exists an integrable dominating function. 

B The vanishing remainder term 

As announced in Remark 15.21 we now study the remainder terms resulting from the expansions 
employed in Section [5.31 Some calculus in normed spaces is needed. 
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B.l Mean Value Theorem 

The key tool is a Mean Value Theorem for Gateaux derivatives. To introduce the central objects 
let Yi be a normed vector space (not necessarily complete) and I2 be a Banach space. Then the 
space L(Yi, Y2) of all bounded, linear functions from Yi to Y2 equipped with the operator norm is 
a Banach space [551 Thm. II. 1.4]. If for a function F : Yi ^ Y2 the strong limit 

, Fix + th) -Fix) , ,, , , 

lim ^ — '- !-^ =: Fu\x]h (B.l 

t->o t 

exists, it is called the Gateaux differential at x in direction h. The Gateaux differential is unique, 
if it exists. We say F is Gateaux differentiate at x if the Gateux differential at x exists in direction 
of all h G Yi and F is Gateaux differentiate , if it is Gateaux differentiable at all x G Yi. If i^ is 
Gateaux differentiable, then we call F^ the Gateaux derivative of F if the map h i-~> Fu[x]h is in 
L(Yi, I2) for every x G Yi. It is not necessarily the case that a Gateaux differentiable F possesses a 
Gateaux derivative. Finally, we say that F is continuously Gateux differentiable when the Gateux 
derivative exists and depends continuously on x, i.e., x 1— >■ F„[x] is a continuous map from Yi into 
the Banach space L{Yi,Y2). Then for continuously Gateaux differentiable F the following Mean 
Value Theorem holds. 

Proposition B.l. Let F be continuously Gateux differentiable, then for every x,h g Yi, 

F{x + h)-F(x) = n F„[2; + C/i]dc)/i, (B.2) 

where the integral in the right hand side is a Bochner integral in the Banach space L{Yi,Y2). 

Proof. We set f{t) := F{x + th) which is Frechet differentiable in t with derivative f'\t\ = Fu\x + 
th]h £ Y2 = i((0, 1), Y2) depending continuously on t. The Fundamental Theorem of Calculus for 
Frechet derivatives implies that the left hand side of (jB.2[) equals 



F{x + h)-F{x) = /i(l)-/i(0) = / F4u + Ch]hdC. 

Jo 

The map ^ H> Fu[x+(^h] is continuous from [0, 1] to L{Yi, Y2) due to the continuous differentiability 
of F for given x,h £ Yi. Hence, the integral in the right hand side of (|B.2p exists in the sense of 
Bochner and its properties allow that 

/ F,[x + Ch] dc) h = [ Fu[x + Ch] h dC. 
Jo ' Jo 



Thus the left and right hand sides in (jB.2p are equal and the proposition is proven. D 

B.2 Remainder terms 

We study the remainder term tB,n in niore detail. The vanishing of tA,n and ep^n follows analo- 
gously. Thus the aim is to show that for all T > 

lim / ^{iyf^h{Z\')tBA^)\dt = 0, (B.3) 

"-*°° Jo 

where tB,^ is the term resulting from the expansion 



The Mean Value Theorem in Proposition IB. II yields an expression for ^B.n resulting from an 
expansion of the map B : H x E ^ X* which is partially continuously Gateaux differentiable due 
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to assumption (CJT]). We obtain the equality 



= B{u,p+{z"-p))+(^j S„[u + C(t/"-u),p+(z"-p)]dc)(C/"-t^) 
= B{u,p) + ( / Bu[u + aU''-u),p+{z^-p)] dC){U''^u) 



eL{H,X*) 



+ (/ Bp[u,p+az"-p)\ d(:){z"~p), 



£L(E,X*) 

where for the sake of simphcity we omit the time argument. By definition (j5.8p the remainder 
term satisfies 



eB,n = {B{U",z"{&')) - B{u,p),^)^ ~ {B^[u,p]{U" - u),^)x - {Bp[u,p]{z"{&') ~ p),ij)x 

Hence we get 

X*3?B.n = [J B.4u + aU''-u),p+iz"-p)]dc){U"^u)-B^[u,p]{U"-u) 

+ (y Sp[«,p + C(z" -P)]dc)(z" -p) - Bp[u,p]{z^~p) . 
Therefore, the remainder term is given by 
i^d^h{Z^)eB,n = ^{[J Bu[u + C{U^-u),p+{z"-p)]-Bu[u,p]dc)[V^,{U"~u)yi^)^h{Z") 

+ t{U Bp[u,p + C{z^-p)]-Bp[u,p]dc)[V^{z"~p)yi;)^h{Z^). 



This is a continuous function of the random vector ([/" — u, z" — p, y/a^{W- — u), y/a^{z^ — u)) 
which converges weakly to the random vector (0, 0, Z)|j Due to the Continuous Mapping Theorem 
the error thus converges weakly to the functions evaluated at the limit which gives zero. We are 
interested in the convergence of the expectations hence due to ^TH Appendix, Prop. 2.3] it remains 
to show the uniform integrability of the errors i^/a^h{Z") eB,n- 

Due to the de la Vallee-Poussin Theorem it is sufficient to show that the q-ih moments are 
uniformly bounded for some q > 1 which we infer from the Lipschitz continuity of the partial 
derivatives in the following. We choose some q € (0, 2). First we obtain the estimate 

E"||»v^/i(Z")eB,„||^ < E"||V^?B,„||i-. IIVII^ 



®The last two components being Z" converge weakly jointly to Z and the first two components converge due to 
Theorem 13.11 (law of large numbers) weakly to the constant 0. Hence joint weak convergence of the whole vector 
holds. 
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and further estimating the stochastic term in the right hand side yields 

E"||V^eB,„|l^. 



< 2E" 
+ 2E' 

< 2E" 



Bu[u + aU''-u),p+iz"-p)]- Bu [u, p] dC 



Bp[u,p + C{z''-p)]-Bp[u,p]dC 



L{£'.X') 



L{X',X') 
(V^||z"-p||£.) 



{y^\\U"^u\\x.) 



\Bu[u + CiU''^u),p+{z''-p)]-B4u,p]\ 



+ 2E" 
< 2L2 E" 
+ 2L^ E" 



Bp[u,p + Ciz"" ~ p)] - Bp[u,p]\\^^^,^^,^d(:{Va^\\z]} ~ p\\£.) 
{\\U" - u\\x + \\z'^ - p\\e){V^^W - u\\j,,f' 
||^"-p|U(V^IN"-p||£.) 



Here L > 1 is a common Lipschitz constant on the derivatives i3„ and Bp. Next using Young's 
inequahty we obtain 



IE"||V^?B,„||?f. 



< CE" 



die/" - u\\x + 11;^" - pWe) 2-9 + f (V^IIC/" - u\\x>) 1 



"kzRiWTTn 



CE" 



¥ll^"-Plli"'+f(V^II^"-PlU*)' 



2q 2q 

< C(E"\\U'' - uf^-^ +E\\z" - pWI-'^ + a,,E^U" - ufx, + a„E"\\z^ - p\\l 



where 23- i 2 for gf — > 1. Moreover the terms a„||L/" — u||^. and a„ ||z" — p|||, arc uniformly 

bounded due to Lemma BTT] and the terms E||t/" — u|| j^T' and E||z"— pj|^^'' are uniformly bounded 
by assumption (C|5]) for some q £ (1, 2). Thus we obtain 

lim E"[iV^/i(Zf)eB,„(t)] =0 Vt > 



and the just derived uniform boundedness over compact intervals implies (jB.Sp due to the Domi- 
nated Convergence Theorem. 
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